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FOREWORD 

This  report  was  prepared  by  the  National  Engineering  Science 
Company,  (NESCO),  Pasadena,  California,  on  Air  Force  Contract 
No.  AF  33(657)-10794,  "Advanced  Air-Breathing  Engine  Study." 

The  work  was  administered  under  the  direction  of  the  Research  and 
Technology  Division,  Wright -Patter son  Air  Force  Base,  Ohio. 

Mr.  Ruben  Canny  was  Project  Engineer  for  the  Research  and 
Technology  Division. 
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Dr.  Raymond  Kushida  was  the  program  manager  for  NESCO. 
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This  report  is  the  Final  Report  on  Contract  AF  33(657)-10794. 
The  contractor's  report  number  is  SN/129. 

This  report  is  classified  "Confidential"  because  it  contains 
data  on  cycle  description,  engine  and  component  performance,  and 
engine  and  component  operating  conditions  pertaining  to  advanced 
air-breathing  engines. 
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ABSTRACT 

This  report  summarizes  analytical  studies  on  supersonic  com¬ 
bustion  ramjet  engines.  The  performance  capabilities  of  components, 
such  as  inlet,  combustor,  and  exit  nozzle,  are  studied  from  theoretical 
considerations  and  their  influence  on  the  engine  performance  is  exam¬ 
ined.  Turbulent  mixing  in  combustor,  nozzle  flows  with  nonuniform 
entrance,  incomplete  recombination  within  the  nozzle,  and  the  influence 
of  nozzle  contour  are  included. 

This  abstract  is  unclassified. 

This  Technical  Document  has  been  reviewed  and  approved 
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1.0  INTRODUCTION 

The  purpose  of  the  Advanced  Air-Breathing  Engine  Study  is  to 
conduct  analytical  inve  stigations  of  means  to  optimize  the  propulsive 
potential  of  the  hydrogen  fueled  supersonic  combustion  ramjet  engine. 
This  is  the  Final  Report  on  this  study. 

This  report  has  been  organized  so  that  each  section  is  essentially 
complete  within  itself.  References  and  figures  pertaining  to  a  particular 
section  will  be  found  at  the  end  of  each  section. 


Manuscript  released  by  the  author  4  August  1964  for  publication  as 
an  RTD  Technical  Documentary  Report. 
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2. 0  SUMMARY 

The  study  of  component  performance  and  their  effect  on  the 
over -all  thrust  performance  of  the  supersonic  combustion  ramjet 
engine  is  clarified  by  explicit  use  of  the  "impulse  function"as  a  variable 
rather  than  pressure  velocity,  enthalpy,  or  temperature.  Consequently, 
the  nonuniformity  at  the  combustor  entrance  due  to  realistic  inlet  con¬ 
figurations  is  shown  to  be  of  secondary  importance  as  far  as  its 
effect  on  the  over -all  performance  is  concerned.  It  is  believed  that 
the  use  of  the  impulse  functions  will  lead  to  significant  simplification 
in  the  application  of  theoretical  computations  and  interpretation  of  experi¬ 
mental  data  to  the  design  and  performance  evaluation  of  engines. 

The  experimental  data  on  penetration  of  fuel  injection  into  a  super¬ 
sonic  stream  could  not  be  correlated  with  theory  based  on  subsonic  ex¬ 
perience.  Supersonic  theories  which  do  correlate  the  experimental  data 
call  for  much  smaller  penetrations. 

The  integral  equations  of  motion  for  turbulent  compressible  mix¬ 
ing  with  equilibrium  reaction  were  used  for  the  study  of  Scramjet  com¬ 
bustion.  Mixing  in  variable  area  ducts  is  studied.  Release  of  heat  due 
to  combustion  shortened  the  flame  length.  Constant  area  ducts  required  a 
shorter  distance  for  the  same  amount  of  mixing  compared  to  divergent 
ducts.  Simulated  shock  losses  on  the  air  side  of  the  mixing  region  caused 
the  turbulent  shear  to  decrease,  leading  to  a  longer  combustor. 

The  thrust  delivered  by  the  divergent  portion  of  the  exit  noszle 
is  studied  with  variation  of  nozzle  configuration  parameters.  For  a 
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nozzle  of  given  exit  area  the  influence  of  wall  friction  leads  to  an 
optimum  length  beyond  which  the  performance  is  reduced.  The  in¬ 
fluence  of  nonuniform  flow  conditions  at  the  nozzle  entrance  and 
incomplete  recombination  within  the  nozzle  are  examined  for  a 
specific  configuration.  It  is  shown  that  these  losses,  even  though 
dependent  upon  the  flow  field  obtained  in  the  nozzle,  can  be  estimated 
from  a  quasi  one -dimensional  analysis. 
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3.  0  ON  SCRAMJET  CYCLE  ANALYSIS 
3.  1  INTRODUCTION 

Parametric  analysis  of  the  performance  potential  of  super¬ 
sonic  combustion  ramjets  (Scramjets)  by  several  groups  (Refs.  1, 

2,  &  3)4  demonstrate  the  sensitivity  of  the  net  thrust  on  component 
efficiencies.  The  strong  interdependence  of  the  component  efficien¬ 
cies  and  the  over-all  thrust  prevents  a  ready  untangling  of  the  relative 
importance  of  the  various  physical  loss  phenomena.  In  the  present 
analysis,  the  engine  cycle  is  analyzed  using  the  impulse  function. 
Concise  and  fundamental  insight  into  the  loss  processes  in  the  Scramjet 
is  thereby  obtained. 

3.2  LIST  OF  SYMBOLS 

A  Area 

Cq  Diffuser  Internal  Drag  Coefficient 

C-,.  Impulse  function  Efficiency 

Cy  Velocity  Coefficient 

D  Drag 

D  Hydraulic  Diameter 

F  Force,  Thrust 

f  Fuel-to- Air  (Masq)  Ratio 

f  Chemical  Heat  Release  Fuel- Air  Ratio 

*  Subject  references  appear  in  Section  3.  11,  page  37. 
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LIST  OF  SYMBOLS  (continued) 


Units 


sf 


Combustor  Friction  Factor 
Enthalpy 

Fuel  Specific  Impulse 
Specific  Impulse  Function  ft/sec  or  lbf- sec/lb 

Kinetic  Energy  Conversion  Efficiency 
Length 


Btu/lb 

lbf  -sec/lb 
i  m 


J 

K 
L 
M 
& 

P 

q 

s 
S 
V 
or 

r 

*lc 

??KE 

nr 

A°f 

Subscripts 
C  Combustor 

D  Diffuser 

f  Fuel 


m 


Mach  Number 
Mass  Flow 
Pressure 

Dynamic  Pressure 
Entropy 
Surface  Area 
Velocity 

Direction  Cosine  Angle 

Ratio  of  Specific  Heats 

Combustion  Efficiency 

Kinetic  Energy  Efficiency 

Static  Pressure  Recovery  Efficiency 

Fuel  Lower  Heating  Value 


ft 


lb  /sec 
m 

lbf  /ft2 
lbf  /ft2 
Btu/lb -°R 
ft2 

ft/sec 

rad 


Btu/lb 


m 
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Subscripts  (continued) 

1  Ideal.  Isentropic 

int  Internal 

N  Exhaust  Nozzle 

t  Total,  Stagnation 

0  Free  Stream,  Ambient 

2  Diffuser  Exit  (Combustor  Entrance) 

3  Fuel  Injection  Station 

5  Combustor  Exit  (Exhaust  Nozzle  Entrance) 

6  Exhaust  Nozzle  Exit 

3.  3  GENERAL  DISCUSSION 

One  of  the  primary  factors  which  distinguish  Scramjet  cycle 
analyses  from  the  conventional  subsonic  combustion  ramjet  cycle 
analyses  is  the  use  of  a  momentum  balance  around  the  combustor 
in  order  to  obtain  the  combustor  exit  velocity  and  pressure.  In  con¬ 
ventional  ramjets  the  combustor  exit  velocity  is  determined  by  the 
exhaust  nozzle  sonic  throat  condition,  hence  is  primarily  a  matter 
of  energy  balance.  The  pressure  at  the  sonic  throat  on  the  other 
hand  is  determined  by  the  inlet  pressure  ratio.  The  coupling  between 
pressure  ratio  and  throat  velocity  is  relatively  weak,  so  that  for  a 
given  combustor,  the  thrust  can  be  directly  related  to  the  inlet  per¬ 
formance  alone.  In  a  Scramjet  however,  the  combustor  exit  velocity 
and  pressure  are  strongly  dependent  on  each  other  c.nd  the  role  of 
inlet  performance  is  not  clear-cut.  Another  way  to  look  upon  the 
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inlet  is  primarily  as  a  means  of  modifying  the  stream  thrust  rather 
than  as  a  thermodynamic  process.  The  use  of  internal  drag  coef¬ 
ficients  rather  than  entropy  rise  as  a  measure  of  losses  suggests 
itself.  Evaluation  of  gas  properties  can  then  be  avoided  in  approxi¬ 
mate  analyses  at  least  to  the  nozzle  expansion  process,  since  the 
drag  coefficients  already  reflect  real  gas  effects. 

The  use  of  various  measures  of  entropy  rise  due  to  irrevers¬ 
ible  losses  in  the  Scramjet  diffuser  has  lead  to  a  number  of  highly 
arbitrary  modifications  of  conventional  efficiency  parameters.  For 
the  sake  of  convenience  in  computation,  the  kinetic  energy  convers¬ 
ion  efficiency,  K^,  and  the  pressure  recovery  efficiency p  have  been 
used  in  place  of  the  more  conventional  kinetic  energy  efficiency  ^ 
and  total  pressure  ratio.  Since  the  entropy  rise  must  be  interpreted 
in  terms  of  thermodynamic  state  variables  such  as  temperature,  pres¬ 
sure,  enthalpy,  and  density,  various  definitions  of  the  "correct  aver¬ 
age  state  properties"  have  been  advanced  for  situations  where  the 
flow  is  highly  nonuniform.  Entropy  changes  occur  due  to  both  vis¬ 
cous  dissipation  (in  boundary  layers  and  in  shock  waves)  and  heat 
losses,  while  stream  thrust  changes  occur  only  due  to  forces  acting 
on  the  control  surface  (body  forces  such  as  bouyancy  are  negligible). 
Arbitrary  definitions  of  efficiency  parameters  are  introduced  to 
insure,  in  an  approximate  manner,  that  the  parameter  applies  to  "an 
equivalent  adiabatic  process".  Unfortunately,  the  interpretation  of 
detailed  component  investigations  in  terms  of  these  entropy-rise 
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efficiencies  are  not  always  consistent  with  the  usage  in*cycle  analyses, 
and  the  performance  figures  obtained  may  be  quite  far  removed  from 
what  actually  is  implied  in  experiments,  or  component  analyses. 

In  the  exhaust  nozzle,  the  large  part  of  the  flow  is  adiabatic 
and  generally  shock  free.  Irreversible  processes  are  either  due  to 
boundary  layers,  which  are  restricted  to  a  relatively  small  part  of 
the  flow,  or  due  to  chemical  nonequilibrium  which  is  best  understood 
from  a  thermodynamic  viewpoint.  Hence,  thermodynamic  measures 
*  of  exhaust  nozzle  losses  are  convenient.  However,  appropriate 
definitions  of  mean  flow  quantities  are  still  troublesome. 

In  the  simplified  model  of  the  engine,  a  fixed  geometry  engine 
with  a  constant  area  combustor  and  an  exhaust  nozzle  exit  area  equal 
to  the  frontal  area  of  the  engine  is  taken  as  a  reference.  In  the  sense 
that  most  engine  designs  will  have  a  fixed  geometry,  themodel  is  some¬ 
what  more  realistic  than  rubber  engines  where  the  pressure  or  en¬ 
thalpy  ratios  are  fixed.  Let  us  consider  an  idealized  supersonic  com¬ 
bustion  ramjet  as  sketched  in  Figure  3.1.  By  convention,  ramjet  inter¬ 
nal  thrust  is  defined  by  the  net  change  in  momentum  of  the  fluid  pass¬ 
ing  through  the  engine,  less  the  ambient  pressure  area  terms: 

Fist  *  '"tX  ~  (3.  la) 

In  practice,  especially  with  the  exotic  nozzle  configurations,  the  exit 
velocity  V^,  pressure  P^,  and  even  the  area  can  be  difficult  to 
define  in  a  manner  suitable  for  engine  performance  calculations.  It 
is  convenient  therefore  to  define  a  mass  flow  average  specific  impulse 
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ENGINE  STATION  LOCATIONS 

0  FREESTREAM 

2  DIFFUSER  EXIT 

3  FUEL  INJECTION 

5  COMBUSTOR  EXIT 

6  NOZZLE  EXIT 


FIGURE  3.  1 

SCHEMATIC  SC  RAMJET 
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function  which  would  be  a  direct  measure  of  stream  thrust  per  unit 
of  mass  flow.  We  shall  designate  the  specific  impulse  function  by 
the  symbol  j  and  define  it  as  an  integral  quantity  averaged  to  take 
into  account  nonuniform  flow:  . 

J  TvxcJm  +  -ft  f  PG>s*<iS 

/t 

where  Vx  is  the  axial  component  of  velocity  and  OC  is  the  angle  of 
the  normal  to  the  elemental  control  surface  with  respect  to  the  axis. 
In  subsequent  discussion  it  will  be  assumed  that  the  velocity  vector 
is  substantially  aligned  with  the  axis,  and  that  the  projected  elemen- 
talareadA  designates  Cos  OC  dS.  Equation  (3.  la)  is  then  written  in 
the  more  general  form: 

fist*  "Uj*  ~  +  P*(A*“M 


where  A^  is  defined  as  the  projected  area  of  the  nozzle  exit  in  a 
plane  normal  to  the  axis,  and  Aq  is  the  capture  area  at  free  stream 
conditions  of  the  air  which  passes  through  the  engine.  The  fuel  to 
air  ratio  will  be  symbolized  by  f,  so  the  net  fuel  specific  impulse 
would  be 


We  shall  modify  the  Eq.  (3.  3)  so  that  ratios  of  the  j  functions  across 
the  inlet,  combustor,  and  exhaust  nozzle  are  called  out: 


Isf 


.it 

Jr  J* 


Ve 


o- 
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The  various  impulse  functions  are  defined  below: 

mass  flow  average  impulse  function  at  the  nozzle 
exit 


j^.  impulse  function  at  the  nozzle  exit  for  isentropic  ex¬ 
pansion  over  the  same  area  ratio. 

jg  mass  flow  average  impulse  function  at  the  combustor 
exit 

mass  flow  average  impulse  function  for  the  air  stream 
at  the  combustor  entrance 

free  stream  impulse  function 


The  breakdown  of  the  impulse  function  in  this  manner  enables 
one  to  evaluate  the  contribution  of  each  component  of  the  engine  to 
the  development  of  the  net  thrust.  In  the  particular  case  where  the 
free  stream  capture  to  exhaust  nozzle  area  ratio  is  1.0,  then  the 
contribution  of  the  last  term  of  Eq..  (3.4)  is  zero. 

The  net  thrust  of  an  engine  can  be  computed  by  adding  and 
subtracting  the  various  forces  which  are  acting  on  the  captured  air 
stream  and  then  Bumming  up  to  obtain  the  net  gain  or  loss  in  the 
impulse  function.  This  process  is  illustrated  for  flight  velocities 
of  8,  16  and  26  thousand  feet  per  second  using  examples  computed 
using  conventional  efficiency  parameters.  The  results  are  given 
in  Figure.  3.  2.  The  impulse  function  based  on  the  fuel  mass  flow  is 
plotted  against  arbitrary  engine  station  points  corresponding  to  Figure 
3.  1.  From  station  0  to  2,  the  drop  in  the  impulse  function  from 
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FIGURE  3.  2 

IMPULSE  FUNCTION  AT  VARIOUS  ENGINE  STATIONS 
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the  free  stream  value  is  due  to  the  inlet  compression  process.  From 
station  2  to  3,  the  fuel  injection  momentum  and  strut  losses  are  in¬ 
cluded.  The  change  in  the  value  of  the  impulse  function  from  station 
3  to  5  is  attributed  to  friction  alone  since  a  constant  area  combustor 
was  assumed.  The  expansion  in  the  exhaust  nozzle  occurs  between 
stations  5  and  6.  The  difference  between  the  fuel-flow-based  free 
stream  impulse,  jQ/f,  and  the  final  value  at  the  end  of  the  exhaust 
nozzle,  (1  +  f)j^/f,  is  the  net  fuel  specific  impulse.  The  impulses 
were  computed  for  two  values  of  the  nozzle  efficiency  parameter  K^. 
The  results  plotted  in  the  form  of  Figure  3.  2  are  very  revealing. 

One  can  immediately  discern  the  effect  of  combustion  friction  losses, 
and  the  relative  importance  of  the  inlet,  fuel  injector,  combustor, 
and  exhaust  nozzle. 

Despite  the  large  differences  in  the  absolute  magnitude  of  the 
impulse  functions,  the  relative  importance  seems  to  be  the  same  for 
each  component  at  all  speeds. 

3.4  INLET  EFFICIENCY  PARAMETERS 

Scramjet  cycle  analyses  usually  start  with  some  arbitrarily 
defined  inlet  efficiency  based  on  thermodynamics.  The  most  com¬ 
monly  used  ones  (Refs.  1,  2,  3,  and  4)  are  kinetic  energy  convers¬ 
ion  efficiency  Kjy  the  kinetic  energy  efficiency*^  the  pressure 
recovery  efficiency  7^  p,  and  the  total  pressure  ratio  P^/P  q. 

These  efficiencies  can  all  be  related  to  the  entropy  rise  associated 
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with  the  shock  and  friction  losses  in  the  inlet.  In  cases  where  the 

air  can  be  assumed  to  be  in  thermal  equilibrium,  the  relations  are 

best  illustrated  on  a  Mollier  chart,  as  shown  in  Figure  3.  3.  The 

initial  conditions  are  at  (h  ,  s  )  where  the  ambient  pressure  is  P 

o  o  o 

and  total  enthalpy  is  The  state  at  the  diffuser  exit  is  (h^, 

where  the  pressure  is  P2  and  the  total  enthalpy  is  h^.  The  difference 
in  total  enthalpy  is  attributed  to  heat  loss  in  the  inlet  stream.  Vari¬ 
ous  hypothetical  state  conditions  are  introduced  in  order  to  obtain 
reference  states  on  which  to  base  an  efficiency  measure.  These 
various  thermodynamic  efficiencies,  defined  in  terms  of  the  quan¬ 
tities  illustrated  in  Figure  3.3,  are 


ht»-  hoi 
hto  -  ho 


(3.5a) 


Hi  -  hoi, 
Kd  *  h'c -  ho 


(3.5b) 


(3.5c) 


Some  variants  of  these  definitions  are  used,  especially  where  there 
is  a  heat  loss  to  be  considered,  but  those  listed  above  are  the  most 
common. 

3.5  INLET  INTERNAL  DRAG  COEFFICIENT 

Inlet  performance  can  also  be  expressed  in  terms  of  aero¬ 
dynamic  coefficients.  Among  these  are  the  internal  drag  coefficient 

of  the  inlet,  and  the  stream  impulse  ratio  j,/j  for  a  fixed  area  ratio. 

c  o 
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hQ  free  stream  enthalpy 

1*2  diffuser  exit  enthalpy 

A  h  heat  loss  in  diffuser 
h^  free  stream  total  enthalpy 

h(2  diffuser  exit  total  enthalpy 

h^  diffuser  exit  enthalpy  with 

heat  loss  restored 

h'  .  enthalpy  for  isentropic  expan - 
ol  sion  from  h',,  P,  to  P 

c  L  O 

PQ  free  stream  static  pressure 

P2  diffuser  exit  static  pressure 

P2j  diffuser  exit  static  pressure 

for  isentropic  adiabatic  com¬ 
pression 

Pto  free  stream  total  pressure 

P^2  diffuser  exit  total  pressure 

SQ  free  stream  entropy 

S2  diffuser  exit  entropy 

S'?  diffuser  exit  entropy  for  heat 
^  loss  restored  at  diffuser  exit 
pressure 


ENTROPY 

•  Actual  State  Points 

Q  Hypothetical  State  Points 

FIGURE  3. 3 

SCHEMATIC  MOLLIER  CHART 
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Cd 


If  the  drag  force  of  the  inlet  is  measured,  then  an  internal  drag 
coefficient  can  be  defined  as 

Ds-Pa 
*1#  Ao 

where  Aq  is  the  free  stream  capture  area,  Dg  is  the  drag  force  on 
the  internal  wetted  surface  and  D&  is  the  additive  drag  due  to  spill¬ 
age  around  the  cowl.  The  ratio  of  the  stream  impulse  function 
is  then  simply  related  to  Cq  by 


£ 


l- 


m 


I  +  TM5I  z  ■  * 

Either  C  or  the  j,/j  ratio  can  be  used  as  measures  of  inlet  per- 
formance.  Note  that  enters  directly  into  the  combustor  momentum 
balance,  hence  pressure,  velocity,  temperature  and  density  are  not 
required  except  to  aid  in  evaluating  combustor  friction,  wall  pres¬ 
sures,  and  heat  loss.  A  plot  of  j~/j  calculated  by  Eq.  (3.7)  is 
given  in  Figure  3,4  for  various  values  of  C^,  and  Mach  numbers. 

The  influence  of  the  area  ratio  A0/A  is  so  small  that  it  is  indetect- 

6  o 

able  on  this  plot. 

The  relation  between  and  the  kinetic  energy  conversion 
efficiency  Kq  are  given  in  Figures  3.5,  3.6,  3.7,  and  3.  8  for  given 
values  of  capture  area  ratio  Aq/A ^  and  flight  Mach  number.  These 
Kp  values  apply  only  to  cases  of  adiabatic  flow,  and  uniform  com¬ 
bustor  inlet  conditions.  Conversion  to  other  efficiency  measures 
can  be  obtained  from  plots  given  in  Refs.  1  and  3.  For  the  one  case 
reported  so  far  where  Scramjet  inlet  drag  has  been  measured 
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FIGURE  3.4 

INLET  PERFORMANCE-  INLET  IMPULSE  RATIO  VS 
INTERNAL  DRAG  COEFFICIENT 
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FIGURE  3. 5 

INLET  INTERNAL  DRAG  COEFFICIENT  VARIATION  WITH  INLET 
KINETIC  ENERGY  CONVERSION  EFFICIENCY 
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KINETIC  ENERGY  CONVERSION  EFFICIENCY,  K0 


FIGURE  3.7 

INLET  INTERNAL  DRAG  COEFFICIENT  VARIATION  WITH  INLET 
KINETIC  ENERGY  CONVERSION  EFFICIENCY 
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experimentally  (Reft.  2,  11),  the  experimental  Cq  at  a  function  of 
area  ratio  are  given  in  Figure  3.9. 

The  experimentally  measured  internal  drag  coefficient 
implicitly  contains  terms  for  the  drag  due  to  the  inviscid  pressure, 
the  induced  pressure  due  to  boundary  layer  displacement,  and  the 
wall  friction.  Corrections  of  experimental  data  on  model  inlets  to 
full  scale  inlets  can  then  be  made  by  applying  corrections  to  each  of 
these  terms.  In  this  respect,  the  hypersonic  simulation  technique 
used  by  the  General  Electric  Company  (Ref.  1)  is  probably  the  best 
method  of  obtaining  applicable  data  short  of  full  simulation  testing. 
From  hypersonic  scaling  relations  (Ref.  10),  the  contribution  of  the 
inviscid  pressure  field  and  the  induced  pressure  to  the  overall  internal 
drag  coefficient  C^  are  the  same  for  the  model  as  for  the  full  scale 
provided  that  the  similarity  parameters  (Md/L),  M^  ‘Yc*  /  Yr,  h  /h 

W  6 

3 

and  M(t/L)  are  held  constant.  There  is  however  a  correction  for  the 
part  of  the  drag  coefficient  due  to  viscous  shear.  The  viscous  shear 
losses  will  increase  approximately  as  the  Mach  number.  Hence  simu¬ 
lated  inlet  drag  coefficient  data  obtained  at  Mach  16  will  be  about  3  %  to 
16%  lower  than  if  the  inlet  were  tested  at  Mach  25.  This  correction 
is  probably  negligible  in  the  first  approximation  but  should  be  included 
if  the  test  data  on  drag  coefficients  are  sufficiently  reliable.  The 
same  correction  should  be  applied  to  the  measured  pressure  recovery 
coefficient  p  although  the  percent  correction  would  vary  with  the 
actual  value  of  ^  p. 
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The  substitution  o£  a  gas  with  a  different  specific  heat  ratio 
such  as  helium  would  not  maintain  similarity.  The  specific  heat 
ratio  is  a  similarity  parameter  in  itself  (Ref.  10),  hence  theory 
requires  that  this  parameter  be  the  same  to  maintain  correspondence. 

The  General  Applied  Science  Laboratory  has  used  the  prin¬ 
ciples  of  hypersonic  similitude  in  nozzle  flow  (Ref.  3).  In  one  test 
setup,  the  geometry  and  specific  heat  ratio  are  simulated  while  an 
entirely  different  test  is  required  to  maintain  chemical  similarity. 

Heat  loss  to  the  walls  of  the  inlet  would  be  expected  to  have  a 
relatively  small  effect  on  the  drag  coefficient  C^.  However  it  can 
have  a  large  effect  on  the  ave  rage  static  pressure  measured  at  the 
diffuser  exit.  Using  experimentally  determined  static  pressures  at 
the  exit  station  can  lead  to  gross  errors  in  the  reported  efficiency 
parameters  unless  the  heat  losses  are  measured  and  inlet  nonuniform¬ 
ity  taken  into  account. 

Theoretical  performance  estimates  for  Scramjets  along  pro¬ 
jected  trajectorie s  have  been  reported  in  Refs.  1  and  2.  The  inlet 
performance  given  in  these  parametric  studies  has  been  interpreted 
in  terms  of  fuel  specific  impulse  as  a  function  of  the  inlet  internal 
drag  coefficient  as  shown  in  Figure  3.  10.  The  linear  decrease  of 
I  I  with  increase  in  Cjj  may  be  noted.  For  comparison  experimental 
inlet  drag  values  are  about  0.  25  to  0.  33  for  approximately  the  same 
Mach  number  and  range  of  area  ratios.  The  discrepancy  between 
the  required  and  experimental  inlet  drags  are  obvious. 
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3.  6  NONUNIFORM  FLOW  AT  THE  DIFFUSER  EXIT 

A  number  of  studies  purporting  to  include  the  effect  of  non- 
uniform  flow  at  the  diffuser  exit  on  performance  have  been  made 
(Refs.  1,  5,  6).  The  study  in  Ref.  6  considers  only  the  boundary 
layer  displacement  and  momentum  defect  as  a  source  of  flow  dis¬ 
tortion.  As  pointed  out  in  Ref.  4  this  is  probably  the  lower  limit 
of  distortion  effects  on  the  mean  static  pressure  ratio.  The  flow  dis¬ 
tortions  considered  in  Ref.  1  assume  again  that  boundary  layers  are 
the  main  disturbance .  Both  the  analyses  in  Refs.  1  and  6  have  in 
common  the  assumption  that  the  flow  distortions  are  essentially 
carried  through  to  the  ends  of  the  combustor  relatively  undisturbed 
by  the  mixing  and  combustion  of  the  fuel.  Detailed  studies  (Ref.  7), 
however,  indicate  that  the  effects  of  combustor  entrance  boundary 
layer  do  not  persist  very  far  into  the  combustor.  The  type  of  non- 
uniform  flow  obtained  at  the  combustor  exit  is  much  more  a  charac¬ 
teristic  of  the  combustor  itself,  and  not  of  the  initial  distortions. 
This  conclusion  will  be  valid  even  if  the  major  cause,  of  inlet  nonuni¬ 
formity  is  due  to  the  inlet  shock  structure  (Ref.  4).  It  follows  there¬ 
fore  that  to  a  fair  approximation,  there  is  no  need  to  know  the  com¬ 
bustor  entrance  flow  nonuniformities.  Only  the  mass  flow  average 
impulse  j^,  total  enthalpy  h^,  and  the  mass  flow  per  unit  area  need 
be  determined  at  the  combustor  inlet.  In  order  to  make  this  point 
clearer,  consider  an  inlet  which  has  a  fixed  drag  coefficient,  as 
defined  in  Eq.  (3.  6)  and  a  fixed  area  ratio.  Then  will  be  fixed. 
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There  are  a  number  of  different  designs  which  would  give  the  same 
value  of  but  due  to  the  different  configurations,  the  profiles  of 
velocity  would  be  different.  If  these  inlets  were  attached  to  the  same 
combustor,  then  it  can  be  assumed  for  not  too  drastic  distortions 
that  the  flow  coming  out  of  the  combustor  will  have  essentially  the 
same  profile.  The  hypothesis  of  independence  of  the  combustor  exit 
flow  profile  from  the  inlet  flow  profile  should  be  more  valid  with 
highly  efficient  mixing  and  large  fuel  flow  rates.  This  does  not  imply 
that  inlet  flow  distortions  are  not  important  since  there  may  be  local 
interaction  of  the  flow  with  the  fuel  injection  which  can  cause  separa¬ 
tion  and  result  in  incomplete  mixing.  It  does  mean  however  that  with 
fairly  efficient  inlets  matched  to  fairly  efficient  combustors,  changes 
in  the  inlet  displacement  thickness  or  similar  nonuniformity  param¬ 
eters  should  have  negligible  effect  on  the  over-all  engine  performance. 

If  the  viewpoint  of  the  above  argument  is  accepted,  it  be¬ 
comes  obvious  that  the  major  purpose  of  any  measure  of  the  Scram- 
jet  inlet  performance  should  be  that  of  producing  the  correct  value 
of  the  impulse  function  at  the  diffuser  exit.  Most  measures  of  inlet 
diffusion  efficiency  (e.g.  ^  K^,  *7^  p,  or  P^/P^)  are  defined 

using  the  pressure  ratio  and  enthalpy  ratio.  It  has  been  shown  how¬ 
ever  (Ref.  4)  that  of  all  the  mean  properties  one  can  find  at  the  dif¬ 
fuser  exit,  the  pressure  is  the  one  most  sensitive  to  flow  distortion. 
Hence,  the  use  of  a  pressure  ratio  in  the  definition  of  an  efficiency 
must  necessarily  be  accompanied  by  a  specification  of  an  inlet  non¬ 
uniformity  parameter.  Otherwise  the  use  of  this  efficiency  value 
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becomes  subject  to  gross  errors. 

To  illustrate  this  point,  we  use  an  analysis  due  to  Ref.  8  for 
a  Mach  16  inlet  with  an  area  ratio  A^/A^  of  25.  In  Figure  3.  11a 
the  impulse  ratio  j ^ / jQ  is  plotted  as  a  function  of  the  internal  drag 
coefficient  C^.  Figure  3.  lib  is  taken  directly  from  Ref.  8  which 
shows  the  relation  between  the  displacement  thickness  the 

pressure  recovery  coefficient  9^  p,  and  drag  coefficient.  If  we  com¬ 
pare  two  inlets  with  the  same  area  ratio  and  drag  coefficient  but  with 
different  displacement  thickness,  different  values  of  the  pressure 
recovery  p  would  be  computed  (see  dashed  lined  in  Figure  3.11) 
yet  the  engine  performance  of  these  two  would  be  exactly  the  same, 

since  the  same  value  of  j7/j  would  be  obtained. 

Ct  o 

3.7  COMBUSTOR  FRICTION 

The  impulse  at  the  diffuser  exit,  j^,  enters  directly  into  the 
momentum  balance  around  the  combustion  chamber.  The  equation 
for  the  combustion  chamber  momentum  balance  in  terms  of  the  im¬ 
pulse  functions  is 

where  is  the  impulse  of  the  fuel  in  the  axial  direction  evaluated  at 

the  combustor  control  surface;  f  ,  q  ,  L.  and  S  are  the  combustor 

c  c  c  c 

friction  factor,  dynamic  pressure,  length,  and  wetted  perimeter, 
respectively,  and  the  last  term  of  Eq.(3.8)  is  the  wall  pressure  inte¬ 
gral  in  the  combustor.  The  equation  for  energy  conservation,  assum¬ 
ing  that  the  radiative  heat  loss  (as  opposed  to  regenerated  heat  loss) 
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is  known,  is 


ht+  +  f>if  +  2Bbl£122l  _ 


where  is  the  sensible  enthalpy  of  the  fuel  before  regenerative 
cooling,  f'  X°f  is  the  effective  heat  release,  and  Ahc  is  the 
permanent  heat  loss  within  the  combustor.  A  seeming  anamoly 
exists  in  evaluating  at  the  combustor  control  surface  where  the 
fuel  temperature  is  at  its  maximum  and  evaluating  the  energy  bal¬ 
ance  with  the  fuel  at  the  entrance  temperature.  This  is  resolved 
when  one  considers  that  the  increase  in  fuel  temperature  occurs 
only  due  to  convective  cooling  and  must  be  accompanied  by  friction 
losses.  With  complete  regenerative  cooling,  the  increase  in  fuel 
impulse  due  to  combustion  chamber  cooling  generally  will  not 
completely  balance  the  loss  in  combustor  exit  impulse  due  to 
friction. 

From  an  approximate  evaluation  of  the  combustor  friction  in 
a  constant  area  supersonic  duct,  the  wall  friction  term  can  be  put 
proportional  to  the  impulse  function  at  the  combustor  exit 

where  T  is  the  mean  combustor  friction  factor  and  D  is  the  hydraulic 
c 

diameter  of  the  duct.  Thus  Eq.(3.8),  for  a  constant  area  duct,  be¬ 
comes 

ir  = 

0+ Lc/D) 


30 


(3.9) 


(3.  10) 


(3. 11) 


CONFIDENTIAL 


APL  TDR  64-112 


CONFIDENTIAL 


The  ratio  of  the  fuel  impulse  to  the  air  impulse  at  the  com¬ 
bustor  entrance  / j ^  can  be  taken  as  approximately  equal  to  the 
velocity  ratio  if  both  are  sufficiently  supersonic  (M^  >  2,  >  2). 

For  typical  conditions,  the  ratio  j^/j^  ranges  from  2  at  Mach  8  to 
about  1/2  at  Mach  25  (Ref.  1  and  2)  when  the  fuel  temperature  is 
held  at  2000°  F.  However,  if  fuel  temperature  is  allowed  to  vary 
with  combustor  heat  loss,  then  will  increase  approximately  as  the 
square  root  of  the  total  regenerative  heat  loss.  The  net  effect  of 
regenerative  cooling  is  that  the  ratio  j^/j^  expressed  in  Eq.  (3.  11) 
will  decrease  below  the  value  obtained  with  no  regeneration. 

3.  8  COMBUSTOR  EXIT  AND  EXHAUST  EXPANSION 

The  main  purpose  of  the  combustion  duct  is  to  mix  the  fuel 
and  air  and  to  obtain  most  of  the  chemical  reactions  of  combustion. 
On  the  other  hand  the  main  purpose  of  the  exhaust  nozzle  is  to  in¬ 
crease  the  momentum  of  the  flowing  stream  as  much  as  possible. 

The  combustor  losses  in  momentum  discussed  previously  are  a  per¬ 
turbing  factor  on  the  operation  of  the  combustor  while  the  major 
processes  of  the  combustor  are  considered  as  perturbations  in  the 
operation  of  the  exhaust  nozzle.  Hence  it  is  convenient  to  separate 
the  considerations  of  combustor  and  exhaust  nozzle,  even  though  the 
same  physical  and  thermodynamic  processes  are  operative.  The 
dividing  control  surface  is  not  always  clearcut,  but  it  will  be  assum¬ 
ed  that  such  a  division  can  be  made. 

In  most  analyses  of  the  exhaust  nozzle  it  is  assumed  that 
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heat  release  due  to  residual  combustion  is  equal  to  that  o£  an  equiva¬ 
lent  equilibrium  gas  down  to  a  freezing  point,  and  also  that  no  further 
mixing  occurs  within  the  nozzle.  The  combustion  is,  therefore, 
charged  with  all  losses  in  heat  release  due  to  incomplete  reaction  and 
mixing.  If  the  combustor  exit  velocity,  enthalpy,  pressure,  and  mix¬ 
ture  ratio  were  known  for  each  stream  tube,  it  is  a  fairly  straight¬ 
forward,  though  tedious,  procedure  to  use  a  rotational  method  of 
characteristics  to  evaluate  the  inviscid  flow  field  of  the  nozzle.  However, 
these  procedures  are  restricted  to  carefully  specified  combustor- 
nozzle  combinations.  Therefore,  to  obtain  more  general  analyses, 
a  means  of  obtaining  useful  efficiency  parameters  is  desirable.  It 
has  often  been  stated  (Refs.  1,  2,  5,  6)  that  it  is  not  possible  to 
choose  mean  values  of  temperature,  pressure,  density,  and  velocity 
to  satisfy  simultaneously  the  requirements  of  mass,  momentum, 
area  and  energy  conservation.  However,  the  viewpoint  taken  at 
NESCO  is  that  there  are  mean  values  which  satisfy  these  require¬ 
ments  provided  that  one  uses  properly  defined  nonuniformity  coef¬ 
ficients  in  the  manner  of  Crocco  (Ref.  12).  The  "ideal"  nozzle 
process  is  defined  in  this  case  as  the  one -dimensional,  isentropic 
expansion  of  a  uniform,  equilibrium  combustion  gas  having  the  same 
combustor  exit  impulse,  fuel-to-air  ratio,  mass  flow  per  unit  area, 
and  total  enthalpy  (after  reaction  and  mixing  heat  losses  are  sub¬ 
tracted)  over  the  given  area  ratio  of  the  nozzle.  The  equivalent, 
real  expansion  process  would  then  have  losses  due  to  nozzle  friction, 
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heat  loss,  nonequilibrium  reactions,  residual  mixing,  contouring, 
and  nonuniform  entrance  conditions.  The  loss  in  performance  is 
then  usually  evaluated  as  a  velocity  coefficient,  kinetic  energy  con¬ 
version  efficiency,  gross  thrust  coefficient,  or  an  impulse  efficiency 
ratio . 

The  velocity  coefficient,  Cv,  is  defined  as 


The  square  of  the  velocity  coefficient  is  the  kinetic  energy  efficiency 
for  the  nozzle.  The  kinetic  energy  conversion  efficiency,  K^,  is 
defined  as 

-  v«.*  -  Vy  =  -  kk 

Kn  V* t  -  V*  h*-  ■  ■ 

These  kinetic  energy  parameters  suffer  the  same  defect  as  the  ana- 
logeous  ones  in  the  inlet,  namely  that  mean  velocities,  pressures, 
and  enthalpy s  at  the  combustor  exit  and  nozzle  exit  must  be  care¬ 
fully  defined  in  order  to  compute  the  corresponding  impulse.  More 
direct  measures  of  the  efficiency  are  the  gross  thrust  ratio 

Ft_  mu  ~  Po  Al 

Ffci.  mujfcc-  p.  At 

or  the  related  impulse  function  efficiency 


These  two  coefficients  can  be  used  directly  in  the  gross  thrust 
equation  (Eq.  3.  lb)  or  the  impulse  equation  (Eq.  3.4).  The  use  of 
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the  impulae  function  efficiency  is  to  be  preferred  because  inter¬ 
nal  nozzle  losses  do  not  depend  upon  the  ambient  pressure.  It  has 
been  shown  in  Ref.  1  that  the  gross  thrust  ratio  F^/F^.  is  equal  to 
the  velocity  coefficient  Cy.  The  relation  between  Cy  and  has 
been  discussed  in  Ref.  3.  The  relation  between  Cgj  and  Cy  is  pre¬ 
sented  in  Figure  3.  12.  Objections  have  been  raised  to  the  use  of  Cy, 

Cm  *  and  F//F,.  because  the  numerical  values  of  these  coefficients 
Ej  6  6i 

are  very  close  to  unity.  Hence  3  to  4  significant  figures  must  be 
carried  in  computations.  In  the  case  of  Cy  such  objections  are  valid 
since  the  alternative  efficiency  measure  is  more  sensitive  to 
losses.  However  in  the  case  of  Cjy  and  F^/F^.,  it  is  felt  that  the 
slight  disadvantage  of  the  number  of  significant  figures  used  in  com¬ 
putations  is  more  than  outweighed  by  the  advantage  that  such  figures 
of  merit  enter  directly  into  the  nozzle  flow  computations,  without 
requiring  further  evaluation  of  nebulously  defined  mean  exhaust  veloc¬ 
ity  and  pressure.  Gas  properties  enter  into  performance  analyses 
only  in  evaluation  of  well  defined  one -dimensional  uniform  isentropic 
expansion.  The  reduction  of  experimental  data  on  nozzle  flow  be¬ 
comes  considerably  easier  in  this  way. 


3.9  CONCLUSIONS 

1.  The  study  of  the  Scramjet  engine  cycle  is  very  much  hampered 
by  a  cumbersome  means  of  reporting  component  performance 
in  terms  of  "an  equivalent  thermodynamic  process  efficiency", 
especially  in  inlet  investigations.  The  use  of  drag  pressure, 
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and  heat  loss  coefficients  would  be  preferable. 

2.  The  explicit  use  of  the  impulse  function  in  the  momentum 
equation  leads  to  considerable  clarification  of  the  role  of 
various  losses  occurring  in  the  Scramjet  cycle.  In  particular, 
the  effect  of  friction  is  more  easily  understood  and  general¬ 
ized. 

3.  Figures  of  merit  for  inlet  and  nozzles  based  on  impulse  func¬ 
tions  such  as  the  internal  drag  coefficient  of  inlets  are  easier 
to  obtain,  use  and  interpret  than  ones  based  on  enthalpy  or 
velocities. 

4.  The  adoption  of  internal  drag  coefficients  as  a  measure  of 
inlet  performance  will  reduce  the  present  emphasis  on  the 
need  to  evaluate  nonuniform  profiles  at  the  combustor  inlet, 
as  far  as  cycle  analysis  is  concerned.  The  need  to  evaluate 
flow  distortion  for  its  effect  on  fuel  mixing  and  combustion 
still  exists. 

5.  Experimental  testing  of  inlet  models  using  the  principle  of 
hypersonic  similitude  will  lead  to  drag  coefficient  figures 
differing  from  scaled  inlets  only  in  a  relatively  minor  cor¬ 
rection  of  the  friction  drag.  Such  testing  should  be  done  with 
simulation  of  static  enthalpy  levels  if  at  all  possible. 

6.  Inlet  testing  with  gases  other  than  air  is  useful  for  study  of 
separation,  transients,  and  transition  and  for  proof  of 
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theoretical  methods.  Direct  application  of  measured  drags  and 
pressures  obtained  during  such  tests  are  not  consistent  with 
similitude  theory. 

3.  10  RECOMMENDATIONS 

1.  Scramjet  investigation  of  inlets  and  exhaust  nozzle  should  report 
their  results  in  terms  of  drag  coefficients  or  equivalent  measures 
of  impulse  change.  Reduction  of  data  to  kinetic  energy  efficiency 
or  related  thermodynamic  measures  should  be  discouraged. 

2.  Explicit  use  of  impulse  functions  in  cycle  analyses  would  lead 
to  better  understanding  of  performance  sensitivity.  Since  all 
existing  Scramjet  cycle  analyses  are  based  on  the  implicit 
use  of  the  impulse  function  concept,  such  a  change  would  lead 
to  simpler  and  more  straightforward  computer  programs. 

3.  Experimental  results  should  be  reported  or  interpreted  as 
much  as  possible  in  terms  of  forces  acting  on  control  surfaces 
rather  than  in  terms  of  the  "equivalent"  process  which  lead  to 
the  same  surface  forces. 
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4.  0  FUEL  INJECTION  AND  COMBUSTOR  ANALYSIS 

4.  1  INTRODUCTION 

To  obtain  good  performance  in  the  Sc  ramjet  engine,  it  is  nec¬ 
essary  that  the  fuel  mix  and  react  with  captured  air  stream.  The  con¬ 
tents  of  this  section  briefly  summarize  the  conclusions  and  opinions 
formed  through  analytic  investigation  regarding  fuel  injection  and 
mixing  in  Scramjet  combustors. 

4.  2  LIST  OF  SYMBOLS 


E 

f'A 

c 

T  J\ 

c 

H 

j5 

m/A 


s* 

m 

*?R 

Ac 


Eddy  Diffusivity 

Standard  Heat  Release  at  the  Local  Fuel- Air  Ratio,  f 

Standard  Heat  Release  at  the  Over-all  Fuel- Air  Ratio,  T 

Static  Enthalpy,  Absolute  Basis 

Impulse  Function  at  Combustor  Exit 

Mass  Flow  Per  Unit  Area 

Static  Pressure 

Turbulent  Reynolds  Number 

Radial  Coordinate 

Velocity  of  External  Stream 

Displacement  Thickness 

Over-all  Combustion  Efficiency 

Mixing  Combustion  Efficiency 

Reaction  Combustion  Efficiency 

Standard  Heat  of  Combustion  Per  Unit  Mass  of  Fuel 
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LIST  OF  SYMBOLS  (continued) 

JX  Fraction  Mass  Flow 
(P  Density 

4.  3  ON  DEFINITIONS  OF  COMBUSTION  AND  MIXING  PERFORMANCE 

CRITERION 

Many  recent  detailed  investigations  of  turbulent  compressible 
mixing  (Refs.  1,  2,  3)P*and  combustion  reaction  analyses  (Refs.  4,  5, 

6)  show  evidence  of  considerable  insight,  ingenuity,  and  labor  in  their 
solutions.  However  these  results  are  usually  not  in  a  form  to  be 
applied  directly  in  most  performance  cycle  analyses,  assuming  com¬ 
parable  boundary  conditions,  without  considerable  labor  on  the  user's 
part.  This  is  due,  in  part,  to  the  lack  of  adequate  standards  on  which 
to  report  results  and  to  evaluate  them  for  more  general  correlations. 

In  this  section  we  shall  attempt  to  describe  the  essential  features  of 
parameters  which  would  best  relate  theoretical  and  experimental 
combustor  studies  to  those  used  in  parametric  performance  analyses. 

All  engine  cycle  computations  utilize  an  efficiency  defined  as 
an  equivalent  heat  loss  from  an  equilibrium  gas,  based  on  the  heating 
value  of  the  fuel  and  the  fuel- air  ratio.  Let  us  consider  the  case  where 
the  composition,  temperature,  pressure,  and  density  have  been  com¬ 
puted  by  a  detailed  chemical  reaction  program  for  a  given  (uniform) 
entrance  condition  to  some  point  where  the  reaction  is  not  yet  in 
equilibrium.  Then  the  same  dynamic  situation  should  be  reproduced 
if  the  comparable  equilibrium  gas  had  the  same  pressure  and  density, 

♦Subject  references  appear  in  Section  4.  10,  page  55. 
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but  with  a  total  enthalpy  decreased  by  a  heat  loss  equal  to  {'-•WfV 
Thus,  we  can  define  the  combustion  reaction  efficiency,  for  an 

arbitrary  reacting  gas  at  some  instant  to  be 

/j,  |  HfP,g,  wow-aQ-  H  (ftft 

‘  r  Ac 

where  the  H  are  computed  on  an  absolute  basis.  Here  f'  is  the  fuel- 
air  ratio  for  fuel-lean  mixtures  and  is  the  stoichometric  ratio  for  fuel- 
rich  mixtures  while  \c  is  the  standard  heating  value  per  unit  mass  of 


(4.1) 


In  the  case  where  there  is  simultaneously  incomplete  mixing 
and  incomplete  reaction,  the  balance  on  the  over-all  total  enthalpy 
across  the  duct  section  and  the  heat  of  reaction  relation  at  the  refer¬ 
ence  temperature  for  ^  leads  to  the  definition  of  the  over- all  heat 
release  efficiency 


T^c  *  ^  00  d/L 


(4.2) 


where  JU  is  the  fraction  mass  flow.  Here  V  Ac  is  the  reference  heat 
release  if  the  mixture  were  uniform  at  the  over-all  mixture  ratio T. 

The  Ac  inside  the  integral  and  in  the  denominator  are  the  same,  hence 
can  be  cancelled.  Note  that  if  the  mixture  is  everywhere  in  local 
chemical  equilibrium,  then  ^  =  1  and  we  have  the  combustion  mix¬ 
ing  efficiency,  ^  given  by 


(4.3) 
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This  relation  states  that  for  the  case  of  instantaneous  local  equilibrium, 
the  combustion  efficiency^  is  equal  to  the  mixing  efficiency  7\  m>  and 
is  mixing  limited.  If  the  mixture  is  everywhere  less  than  stoichiometric 
and  the  over-all  T  is  less  than  stoichiometric,  or  alternatively,  if  the 
mixture  is  everywhere  greater  than  stoichiometric  and  F  is  greater  than 
stoichiometric  it  follows  that 

7  “X f  <4- 

and  the  mixing  efficiency  Tp  m  is  unity.  Equation  (4.  4)  implies  that  the 
combustion  efficiency  is  now  reaction  limited.  An  interesting  mathemat¬ 
ical  theorem  derived  in  Ref.  7  (page  309)  may  be  applied  to  Eq.  (4.  2) 
for  the  case  of  hydrogen-air  mixing.  Generally,  it  is  the  capture  of 
air  into  the  turbulent  mixing  zone  which  limits  the  uniformization  of 
the  mixture.  Hence,  we  can  say  roughly  that  at  a  given  duct  length 
the  gases  in  stream  tubes  close  to  the  fuel  jet  centerline  have  had  a 
longer  time  to  react  than  those  in  stream  tubes  close  to  the  jet  bound¬ 
ary.  The  reaction  efficiency  should  be  higher  for  highax  fuel-air 
ratios  f.  This  would  make  R  covariant  with  f'  in  the  sense  that 
Crocco  uses  the  term  and  it  follows  using  his  theorem  that 


%« ?  t  ( £  (4' ! 

In  other  words,  the  over-all  combustion  efficiency  ^  c  is  higher  than 
the  product  of  the  mean  reaction  efficiency  and  the  mixing  efficiency. 

Thus,  conservative  estimates  of  ^  are  obtained  by  separately  com¬ 
puting  reaction  and  mixing  efficiency,  and  multiplying  the  two  efficiencies. 
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Notice  that  the  reaction  and  mixing  efficiencies  are  defined  in 
such  a  manner  that  they  are  equivalent  to  local  heat  losses  in  an  equi¬ 
librium  gas.  There  are  other  losses  associated  with  incomplete  re¬ 
action  and  mixing  which  occur  in  the  subsequent  expansion  in  the  nozzle. 
These  may  be  classed  as  nonuniformity  effects.  These  losses,  how¬ 
ever,  are  not  assigned  to  the  combustor,  but  are  expansion  losses 
properly  assigned  to  the  exhaust  nozzle,  as  will  be  discussed  in  Sec¬ 
tion  5.  . 


4.  4  ON  NONUNIFORMITY  AT  THE  COMBUSTOR  EXIT 

There  are  losses  in  the  combustor  due  to  friction  and  heat 
transfer  to  the  walls,  in  addition  to  the  combustion  energy  losses  dis¬ 
cussed  above.  The  wall  friction  is  used  in  the  momentum  balance 
along  with  the  pressure-area  term  (cf. ,  Section  3.6,  Eq.  (3.8))  to 
evaluate  the  combustor  exit  impulse  function  jj. .  The  effective  total 

enthalpy,  H  ,  is  evaluated  from  the  input  total  enthalpy  less  heat  loss 
l5 

and  over- all  combustion  efficiency  loss  defined  as  in  Eq.  (4.2).  The 
quantities  j,.,  Ht  ,  and  rh^/A^  are  sufficient  to  characterize  the  ideal 
thrust  developed  in  an  isentropic  expansion.  Note  that  we  choose  to 
define  quantities  which  would  be  constant  in  an  ideal  constant  area 
combustor  duct.  Hence,  these  quantities  are  chosen  as  the  reference 
constants  upon  which  to  base  the  studies  of  losses.  Had  we  chosen  any 
of  the  quantities  such  as  mean  velocity,  kinetic  energy,  pressure, 
density,  or  enthalpy,  the  evaluation  in  an  ideal  constant-area  combus¬ 
tor  would  require  the  use  of  the  heat,  mass,  and  momentum  balance, 
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implicitly  using  the  quantities  j^,  ,  and  rhg/Ag  but  also  requiring 

5 

further  evaluation  and  definitions  of  the  meaning  of  mean  values.  This 
extra  step  in  analysis  leads  to  some  confusion  regarding  nonuniformity 


corrections. 


Once  we  have  accepted  that  the  basis  for  comparison  is  with 
the  same  exit  conditions  in  the  mean  values  of  impulse  function,  total 
enthalpy,  and  area,  then  the  question  as  to  mean  values  becomes  simply 
a  perturbation  analysis.  The  approach  used  at  NESCO  (Ref.  8)  is 
consistent  with  this  viewpoint.  Other  nonuniformity  analyses,  espe¬ 
cially  those  using  displacement  thickness  in  analogy  to  boundary  layers, 
can  be  easily  adapted. 


4.  5  FUEL  INJECTION 

In  a  preliminary  investigation  into  simplified  theories  of  jet 
penetration,  experimental  data  (Ref.  9)  on  the  cross  stream  penetration 
of  a  supersonic  jet  injected  into  a  supersonic  flow  were  compared  with 
the  theories  of  Lindley  (Ref.  9),  Ferrari  (Ref.  10),  and  Abramovitch 
(Ref.  11).  The  first  two  theories  were  for  supersonic  flow  while  the 
last  was  based  on  subsonic  theory  and  data.  In  every  case,  theory 
called  for  penetrations  which  were  far  greater  by  a  factor  of  at  least 
2  than  was  actually  measured.  Physically  unrealistic  drag  coefficients 
of  the  order  of  50  or  more  were  required  to  rationalize  theory  with 
data.  The  basic  erroneous  assumption  which  lies  at  the  heart  of  each 
of  the  above  mentioned  theories  is  that  the  injectant  stream  is  coherent 
and  quasi-uniform  across  its  diameter.  The  actual  phenomena  calls 
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for  an  expansion  fan  on  the  downstream  side  of  the  jet.  The  theory  of 
Zukoski  and  Spaid  (Ref.  12)  does  not  assume  any  mechanism  for  turn¬ 
ing.  The  injected  stream  immediately  turns  to  flow  along  the  wall. 

The  correlation  with  their  own  data  is  good  and  comparison  with 
Marquardt  data  (Ref.  9)  is  fair.  Since  the  Marquardt  work  is  done 
under  simulated  Scramjet  conditions,  the  latter's  empirical  correla¬ 
tion  seems  adequate  for  most  purposes. 

The  data  of  Zukoski  and  Spaid  on  concentration  profiles  behind 
a  normally  injected  stream  suggest  that  a  vortex  filament  is  shed  from 
both  sides.  The  axis  of  the  vortex  is  roughly  parallel  to  the  wall  and 
has  a  rotation  such  that  material  near  the  wall  is  swept  toward  the 
center  of  the  injected  stream.  The  slot  injector  tried  by  General 
Electric  (Ref.  13),  with  the  long  side  aligned  to  the  flow,  should  cause 
an  even  stronger  pair  of  vortex  filaments  to  be  formed.  The  slot  injec¬ 
tor  obtained  better  mixing  than  the  circular  injector  for  the  same 
pressure  and  flow  rate.  These  vortex  filaments  increase  the  amount 
of  initial  mixing,  but  the  influence  also  seems  to  die  down  rather 
rapidly  since  the  data  reported  in  Ref.  5  with  flush  wall  injection  at 
Mach  2.  3  show  little  improvement  in  mixing  with  distance  downstream. 
These  preliminary  experiments  suggest  that  shaped  slot  injectors  can 
be  designed  which  would  create  a  stronger  single  whirl,  rather  than  a 
pair,  hence  leading  to  more  persistent  and  effective  vortex  mixing. 

At  small  angles  of  inclinations  of  the  injected  stream  to  the 
main  stream,  the  vortex  filaments  do  not  appear  to  form  and  pene¬ 
tration  is  better  (Ref.  13).  It  is  conjectured  that  strong  vortex 
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formation  will  be  associated  with  the  oblique  shock  formed  near  the 
base  of  the  injected  jet.  If  this  shock  wave  is  a  strong  oblique  wave, 
rather  than  a  weak  attached  oblique  shock  then  there  seems  a  greater 
tendency  to  form  vortices.  Hence,  the  amount  of  initial  mixing 
induced  by  the  act  of  injection  should  be  quite  different  when  correlated 
with  the  effective  deflection  angle  of  the  forward  face  of  the  jet. 

The  expansion  of  the  injected  fuel  stream,  where  the  fuel  injec¬ 
tor  pressure  is  much  higher  than  the  combustor  entrance  air  static 
pressure,  will  cause  oblique  shock  waves  to  form.  If  the  blockage  is 
sufficiently  large  then  separation  in  the  boundary  layer  ahead  of  the 
injection  point  will  occur.  Such  separation  by  itself  does  not  affect  die 
over-all  combustor  momentum  balance  and  over-all  performance  unless 
it  affects  pressure  on  convergent  or  divergent  walls. 

Downstream  of  the  injection  point,  certain  local  restrictions 
on  allowable  fuel  injection  conditions  may  occur.  They  are  that  the 
over-all  mass,  momentum,  and  energy  balance,  either  before  or 
after  heat  release,  do  not  reach  choking  conditions.  Certain  other 
limiting  conditions  previously  reported  in  Refs.  8  and  14  were  based 
on  the  assumption  of  uniform  pressure  across  the  combustor.  These 
limiting  conditions  were  found  to  be  easily  relieved  if  the  pressures 
are  allowed  to  be  unequal  in  the  unmixed  or  partially  mixed  fuel  and 
air  stream,  hence  are  not  true  limits  on  injector  configurations. 

4.  6  MIXING  STUDIES 

The  problem  of  mixing  of  supersonic  jets  have  been  investigated 
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in  a  series  of  studies  (Refs.  8,  14,  15,  16).  The  starting  point  for  all 
of  these  analyses  were  the  boundary  layer  assumptions.  As  reported 
in  Ref.  16  the  basic  assumption  of  the  cross  stream  pressure  gradients 
being  negligible  is  suspect,  but  there  is  no  simple  method  for  including 
this  in  mixing  analyses.  The  assumption  that  cross  stream  pressure 
variations  have  a  secondary  effect  on  the  local  spreading  angle,  relative 
to  the  local  velocity  vector  seems  a  reasonable  one  (Ref.  8).  It  is  con¬ 
jectured  that  analyses  that  do  include  the  pressure  effects  will  show 
increases  in  the  over-all  mixing  rate,  but  the  effect  would  be  hard. to 
separate  from  injection  effects. 

The  investigations  have  been  centered  mainly  around  the  use  of 
the  integral  method  of  solving  free  mixing  problems.  In  addition  to  the 
usual  equations  for  the  conservation  of  mass,  momentum,  fuel  mass, 
and  energy,  the  integral  equation  formed  from  the  moment  of  fuel  mass 
conservation  was  used  to  determine  the  rate  of  spread  (Ref.  14).  It 
is  necessary  to  assume  profile  distribution  functions  for  velocity, 
total  enthalpy  and  fuel  mass  fraction,  and  that  this  profile  function  is 
similar  throughout  the  length  of  the  jet.  Based  on  supersonic  wall 
boundary  layer  experience,  it  was  assumed  that  the  proper  basis  for 
similarity  is  the  Howarth  radial  coordinate,  and  on  the  basis  of 
Abramovitch's  experience  in  subsonic  mixing  that  Schlichting 1  s  wake 
profile  function  was  applicable  (Ref.  14).  Comparison  with  data  from 
the  Marquardt  Corporation  (Ref.  5)  and  the  General  Electric  Company 
(Ref.  13)  showed  that  the  first  assumption,  while  requiring  some  fur¬ 
ther  checks,  is  substantially  correct,  but  that  the  second  assumption, 
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coupled  with  the  first,  is  wrong.  The  profile  function  for  hydrogen 
mixing  with  air  in  the  absence  of  reaction  is  considerably  closer  to 
being  Gaussian  than  the  polynomial  expression  of  Schlichting.  For¬ 
tunately,  it  turns  out  that  rather  crude  forms  of  the  profile  function 
lead  to  fairly  good  approximations  to  the  spreading  rate.  This  is,  in 
fact,  one  of  the  primary  advantages  of  the  integral  formulation. 

Solutions  based  on  the  vonMises  transformation  (Ref.  1,  2) 
were  not  used  since  the  solutions  corresponding  to  small  external 
velocity  are  not  correct,  hence  could  not  be  compared  to  most  of  the 
subsonic  data  on  jet  spreading. 

Various  forms  of  the  eddy  diffusivity  formulation  were  investi¬ 
gated.  The  basic  form  which  seems  adaptable  to  both  low- speed  wall 
boundary  layers  and  to  mixing  of  fluids  of  different  density  is  that 
based  on  a  displacement  thickness.  It  is  provisionally  defined  as 

*  Si£r-  *  -jq  (<>«.-  *«)*  (. 

Here  R(  is  a  constant  which  in  two-dimensional  subsonic  jets  takes  on 
the  value  62.  The  form  is  not  inconsistent  with  that  proposed  by  Ferri 
(Ref.  17),  but  is  in  the  form  of  an  integral  rather  than  an  over-all 
mean  difference.  It  has  been  shown  (Ref.  14)  that  the  Eq.  (4.6)  reduces 
to  Ferri's  form  for  low  speed  nonreactive  flow,  and  to  Prandtl's 
classic  relation  for  the  exchange  coefficient  for  constant  density  flow. 

Most  mixing  studies  have  assumed  that  longitudinal  pressure 
gradients  are  negligible.  However,  in  constant- area  supersonic  com¬ 
bustion  ducts,  the  pressure  rise  can  be  quite  large,  especially  at 

48 


CONFIDENTIAL 


APL  TDR  64-112 


CONFIDENTIAL 


lower  flight  speeds.  A  theoretical  study  to  include  these  effects  of 
pressure  gradient  in  a  variable  area  duct  are  reported  in  Ref.  18. 

A  comparison  between  theory  and  data  obtained  by  the  Marquardt 
Corporation  (Ref.  9)  in  cold  flow  tests  on  the  mixing  of  hydrogen  with 
nitrogen  are  shown  in  Figure  4.  1.  Mole  fraction  of  hydrogen  at  the 
centerline  is  plotted  as  a  function  of  the  distance  downstream.  Exam¬ 
ination  of  the  figure  shows  that  the  predicted  centerline  distribution  is 
somewhat  flatter  than  the  experimentally  observed  distribution.  Better 
correlation  could  be  obtained  by  using  a  slightly  smaller  value  of  the 
turbulent  Reynolds  number,  R^.  Considering  that  the  Reynolds  number 
was  evaluated  for  subsonic  jets  of  air  in  a  quiescent  atmosphere,  the 
comparison  is  regarded  as  being  reasonably  good. 

It  was  found  that  mixing  is  more  rapid  in  constant  area  ducts  than 
in  constant  pressure  ducts,  as  shown  in  Figure  4.  2.  Here,  half -width 
of  mole  fraction,  mass  fraction,  and  total  width  are  compared  for  two 
different  assumptions  regarding  the  pressure-area  distribution.  Due  to 
the  density  term  appearing  in  the  eddy  viscosity  formula  (Eq.  (4.  6))  the 
increase  in  pressure  with  mixing  increased  the  eddy  diffusivity,  hence 
decreased  the  mixing  distance.  In  another  case  the  effect  of  shock-waves 
was  simulated  by  using  a  total  pressure  drop  in  the  outer  stream.  It 
will  be  noted  that  the  distance  required  for  the  mixing  zone  to  reach  the 
wall  is  increased  with  total  pressure  loss  in  the  outer  stream,  as  shown 
in  Figure  4.  3.  This  would  indicate  a  combustor  length  trade-off  between 
oblique  shock  waves  caused  by  the  fuel  injector,  and  the  shear -controlled 
mixing  rate.  However,  the  increased  mixing  due  to  the  injector, 
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FIGURE  4. 1 

CENTERLINE  FUEL  DISTRIBUTION 
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which  is  not  considered  here,  and  offset  this  loss. 

The  integral  theory  of  free  mixing  analysis,  while  convenient 
and  rapid,  is  not  easily  adapted  to  the  study  of  simultaneous  mixing 
with  finite  rates  of  chemical  reaction.  The  integral  analysis  requires 
that  the  profile  functions  of  the  reacting  species  be  known.  While  a 
general  polynomial  function  can  be  assumed  with  coefficients  to  be 
evaluated  as  part  of  the  solution,  reliable  results  would  require  a  high 
order  polynomial*  The  direct  numerical  solution  of  the  controlling 
partial  differential  equations,  while  time  consuming,  would  probably 
be  easier  to  program  for  computers.  The  validity  of  such  solutions, 
however,  depend  very  strongly  on  the  validity  of  the  eddy  diffusivity 
relation.  It  would  be  of  interest  to  investigate  some  simple  configura¬ 
tions  for  comparison  with  experiments. 

4.  7  CONSTANT  PRESSURE  MIXING  STUDIES 

Fundamental  studies  on  the  possible  transformations  of  the 
equations  of  compressible  turbulent  mixing  with  reaction  to  that  of  an 
analogous  incompressible  mixing  problem  were  reported  in  Ref.  14. 

In  brief,  it  has  been  shown  that  transformation  of  the  partial  differential 
equations  of  motion  are  not  possible  except  with  physically  implausible 
assumptions  regarding  the  shear  stress  term.  However,  it  is  possible 
to  obtain  a  transformation  of  the  integral  form  of  the  equations  of 
motion.  As  a  consequence  it  is  theoretically  possible  to  correlate 
low-speed  diffusion  flame  data  with  supersonic  data.  Proof  of  the 
validity  of  the  transformation  is  obtained  in  the  prediction  of  core 
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length  in  axisymmetric  jets,  as  was  shown  in  Ref.  14. 

The  transformation  theory  was  applied  to  the  problem  of  a 
two-dimensional  constant  pressure  diffusion  flame  with  equilibrium 
gas  properties.  These  results  were  reported  previously  in  Ref.  14. 
An  example  was  computed  with  air  at  a  velocity  of  3500  fps,  hydrogen 
fuel  injected  at  3000  fps,  air  and  hydrogen  static  temperature  of 
1650°R,  and  1000C  *R  respectively,  and  a  static  pressure  of  2  atmos¬ 
pheres.  Exact  equilibrium  gas  properties  were  used  for  this 
computation.  Profiles  of  velocity,  composition,  temperature,  and 
Mach  number  were  obtained.  Due  to  the  length  of  time  required  to 
compute  one  case  (70  minutes  machine  time),  no  further  examples 
were  computed. 

4. 8  CONCLUSIONS 

1.  Neither  theoretical  nor  experimental  data  on  reaction  and 
mixing  were  found  to  be  reporting  data  in  a  manner  consistent 
with  usage  in  cycle  analyses.  A  definition  of  combustion 
efficiency  as  equivalent  heat  loss  in  an  equilibrium  gas  at 

the  same  pressure  and  density  is  proposed. 

2.  Fuel  jet  penetration  into  supersonic  streams  cannot  be 
correlated  with  theories  which  depend  on  the  drag  to  turn 
the  injected  stream.  Theories  which  do  not  depend  on 
drag,  but  on  Prandtl-Meyer  expansion  of  the  injected 
stream  seems  more  successful. 


54 


CONFIDENTIAL 


APL  TDR  64-112 


CONFIDENTIAL 


3.  Integral  methods  of  analyses  for  mixing  in  variable  area 
ducts  have  been  developed.  Constant  area  ducts  exhibit 
more  rapid  mixing  than  divergent,  constant  pressure  ducts. 
Fairly  good  correlation  with  experimental  data  have  been 
obtained. 

4.9  RECOMMENDATIONS 

1.  Further  work  on  correlating  integral  mixing  analyses  with 
experimental  data  in  constant  area  ducts  is  required. 

2.  A  more  general  analyses  to  include  effects  of  wall  boundary 
layer  and  multiple  injection  is  possible  and  should  be  under¬ 
taken. 
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5.  0  EXHAUST  NOZZLE  STUDIES 

5.  1  INTRODUCTION 

The  nozzle  for  a  supersonic  combustion  ramjet  presents  ana¬ 
lytical  problems  that  are  somewhat  different  from  those  encountered  in 
conventional  rocket  nozzles.  Due  to  the  supersonic  velocity  at  the  nozzle 
entrance,  the.  thrust  increment  obtained  from  such  a  nozzle  is  smaller 
than  the  comparable  conventional  nozzle  with  sonic  throat  velocity. 
Velocity  and  enthalpy  nonuniformities  across  the  combustor  exit,  due 
to  incomplete  mixing  of  fuel  and  air,  are  no  longer  suppressed  by  thermal 
choking.  In  addition  to  the  problems  arising  from  the  supersonic  velocity, 
high  combustion  temperatures  are  responsible  for  dissociation  of  the 
chemical  species  in  the  gases  entering  the  nozzle.  The  resulting  chem¬ 
ical  recombination  and  heat  release  in  the  nozzle  makes  the  gasdynamical 
analysis  of  the  flow  more  uncertain.  The  nozzle  design  will  probably 
be  asymmetric  because  of  unorthodox  inlet  designs  and  over -all  engine 
envelope  requirements. 

The  theoretical  studies  conducted  at  NESCO  on  nozzle  perform¬ 
ance  are  presented  in  this  section  and  can  be  summarized  as, 

a)  examination  of  nozzle  contours  of  limited  length 
designed  for  optimum  thrust 

b)  evaluation  of  loss  in  thrust  due  to  nozzle  wall 
friction 

c)  location  of  "sudden-freezing"  line  in  the  nozzle 
and  its  influence  on  nozzle  thrust  performance 
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d)  changes  in  nozzle  thrust  performance  due  to 
nonuniformities  at  nozzle  entrance 

e)  evaluation  of  lift  force  on  the  nozzle  due  to 
asymmetric  contour  or  nonuniform  entrance 
flow  conditions. 

The  theoretical  approach  used  in  nozzle  flow  analysis  and  discus¬ 
sion  of  the  computational  results  are  given  in  the  following  subsections. 

5.  2  LIST  OF  SYMBOLS 

A  Nozzle  Cross-sectional  Area 

dA  Elemental  Area  of  Control  Surface 
C.,.  Exit  Impulse  Efficiency 

Cj^j  Nozzle  Impulse  Coefficient 

h  Enthalpy 

j  Stream  Thrust  Per  Unit  Mass  Flow 

m  Mass  Flow  Rate  Through  Nozzle 

dm  Elemental  Mass  Flow  Through  Elemental  Area  of 
Control  Surface 

P  Pressure 

s  Coordinate  Along  Stream  Line 

V  Velocity 

X  Axial  Coordinate 

Y  Nozzle  Cross-section  Height  Also  Equal  to  Cross-section 
Area  A  in  the  Case  of  Two-dimensional  Nozzles 

6  Isentropic  Exponent 

y  Velocity  Direction 
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Subscripts 

a,  b,c  Average  Flow  Conditions  Across  Each  of  the  Stream  Tubes 
a,  b,  and  c  (see  Figure  5.8) 

i  Reference  Value  Reached  by  One -dimensional  Flow 

5  At  Nozzle  Entrance  (also  equal  to  subscript  "throat") 

6  At  Nozzle  Exit 

5.  3  NOZZLE  PERFORMANCE  PARAMETERS 

The  primary  purpose  of  a  nozzle,  which  is  to  impart  momentum 
to  the  exhaust  gases  as  they  leave  the  nozzle  exit,  can  be  evaluated  by 
considering  the  exit  impulse.  Let  us  denote  by  subscript  6  the  flow 
conditions  along  a  control  surface  at  the  nozzle  exit  and  normal  to  the 
thrust  axis  as  indicated  in  Figure.  5.  1.  Because,  in  general,  the 
nozzle  length  is  limited,  the  flow  conditions  at  the  exit  control  surface 
would  be  nonuniform  and  one  can  write 


Nozzle  exit  stream  thrust  = 


/Ac  /- m 

pdA  +/  V  Cos  ©dm 


(5 


where 


P  =  static  pressure 
dm  =  elemental  mass  flow 
V  =  velocity 
0  =  flow  direction 


across  an  elemental  area  of  dA  of  the  control  surface.  The  integration 
is  taken  along  the  control  surface  across  the  nozzle  exit  area  Ag  and 
encompasses  the  entire  mass  flow  m  through  the  nozzle.  The  above 
Eq.  (5.  1)  gives  us  the  component  of  exit  stream  thrust  in  the  thrust 
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FIGURE  5.  1 

Schematic  of  exit  nozzle  configuration 
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axis  direction. 

Similarly,  one  can  write 


Nozzle  entrance  stream  thrust  = 


fpd/uJvG«edm 


(5.  2) 


where  the  integration  is  taken  along  a  control  surface 
at  the  nozzle  entrance  area  Ag 


Mass  flow  of  exhaust  gases  crossing  the  control  surface  at  nozzle  entrance 
should  be  equal  to  mass  flow  across  the  control  surface  at  nozzle  exit. 


Hence 


(5.3) 


It  is  found  convenient  to  use  stream  thrust  per  unit  mass  flow  as 


a  parameter  denoted  j.  Hence 

J‘  =  i[PdA  +  f^C°s6d^ 

jf.  +£vCoSedmj 


The  influence  of  the  divergent  portion  of  the  exit  nozzle  is  to  increase 
the  stream  thrust  parameter  from  j  g  at  nozzle  entrance  to  j  g  at  nozzle 
exit. 

For  prescribed  nozzle  entrance  conditions  the  value  of  j  g  depends 
upon  the  properties  of  exhaust  gases  and  nozzle  geometry  such  as  expan¬ 
sion  ratio,  length,  and  wall  contour.  The  latter  effect,  i.e.  ,  the 
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influence  of  nozzle  geometry,  can  be  examined  by  considering  nozzle 
impulse  coefficient  expressed  as 


(5.6) 


In  considering  the  above  coefficient  one  should  remember  that  j  g  is  held 
constant.  The  exit  impulse  function  j  ^  can  be  computed  by  analysis  of 
supersonic  flow  in  the  nozzle  utilizing  an  appropriate  constant  value  of 
^  to  represent  expansion  process  of  the  exhaust  gases.  The  effect  of 
nozzle  wall  friction  also  can  be  included  in  the  evaluation  of  j^.  On  the 
other  hand,  the  influence  of  flow  nonuniformities,  frictional  and  non¬ 
equilibrium  losses  can  be  brought  into  focus  by  examining  the  exit 
impulse  efficiency  defined  as 


Cej  = 


■It. 


(5.  7) 


In  the  above  equation,  j  z  is  the  one -dimensional  value  for  equilibrium 

i 

expansion  of  exhaust  gases  from  uniform  conditions  at  nozzle  entrance 
to  the  particular  value  of  At/A,  under  consideration. 

The  above  defined  exit  impulse  efficiency  is  suitable  to  study 
the  performance  of  a  given  nozzle,  whereas  the  nozzle  impulse  coef¬ 
ficient  as  defined  by  Eq.  (5.  6)  takes  into  account,  also,  changes  in 
nozzle  area  ratio. 

However,  as  can  be  seen  in  the  following  subsections,  the  nozzle 
impulse  coefficient  charts  can  be  easily  prepared  from  nozzle  flow 
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analyses  carried  out  at  an  appropriate  constant  value  of  K  .  The  more 
complicated  computations,  including  nonuniform  entrance  conditions  and 
nonequilibrium  losses,  can  be  evaluated  only  in  terms  of  a  prescribed 
nozsle  geometry. 

5.4  PARAMETRIC  STUDY  OF  NOZZLE  IMPULSE  COEFFICIENT 

Depending  upon  the  SCRJ  engine  and  vehicle  configuration,  the 
exit  nozzle  could  be  simple  conical,  axially  symmetric,  asymmetric 
two-dimensional,  or  three -dimensionally  asymmetric.  To  obtain 
representative  values  for  performance  and  its  dependency  on  nozzle 
area  ratio  and  length, two -dimensional  asymmetric  nozzle  contours  are 
examined.  The  lower  wall  of  the  two-dimensional  nozzle  is  considered 
plane  and  parallel  to  the  thrust  axis,  i.  e.  ,  combustion  chamber  center- 
line.  The  upper  wall  of  the  nozzle  between  the  nozzle  entrance  and  exit 
is  contoured  so  as  to  yield  optimum  thrust.  In  the  case  of  two-dimensional 
nozzles,  Rao's  method  (Ref.  1)*  for  contouring  the  wall  for  optimum 
thrust  becomes  simple  in  that  they  form  truncated  portions  of  contours 
designed  for  appropriate  uniform  exit  Mach  numbers.  In  the  case  of 
axisymmetric  nozzles  such  is  not  the  case  and  wall  contour  computations 
have  to  be  carried  out  for  each  individual  area  ratio  and  length.  It  is 
for  this  reason  that  the  performance  evaluation  given  here  is  based  on 
two-dimensional  asymmetric  nozzles. 

Nozzle  thrust  performance  has  been  computed  for  several 

families  of  two-dimensional  optimal  nozzles  for  two  nozzle  entrance 

conditions  encountered  on  a  typical  flight  trajectory.  The  trajectory 
▼Subject  references  appear  in  Section  5. 12,  page  109. 
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points  and  associated  nozzle  entrance  conditions  were  furnished  by  the 
Marquardt  Corporation.  Appropriate  data  are  listed  below  for  the  two 
points.  The  first  is  a  moderately  low  point  along  the  trajectory  and  the 
second  is  near  the  upper  end. 


Free  stream  Conditions 


Speed  (ft /  sec) 

Altitude  (ft) 

Nozzle  Entrance  Conditions 

15,000 

130,000 

26, 000 

166,000 

Mach  number 

3.  31 

4.  02 

Pressure  (atm) 

1.  064 

0.444 

Temperature  (°R) 

5,  640 

7,  360 

Isentropic  exponent 

1.  105 

1.  179 

Optimal  inviscid  nozzle  contours  were  machine  computed  by  the 
method  of  characteristics  using  Rao's  optimizing  technique. 

The  nozzle  flow  field  computations  were  carried  out  utilizing 
the  method  of  characteristics.  The  details  of  the  procedures  employed 
are  given  in  Ref.  2.  Application  of  this  method  to  nonequilibrium  flow 
of  exhaust  gases  is  much  too  complicated  to  be  of  use  in  parametric 
evaluation  of  nozzle  configurations.  Hence  the  exhaust  gases  are  assumed 
to  behave  as  ideal  gas  with  the  appropriate  values  of  ^  listed  above. 

For  similar  reasons,  the  flow  conditions  at  the  nozzle  entrance  section 
are  considered  uniform  for  the  parametric  study.  However,  the  influ¬ 
ences  of  nonuniformities  at  entrance,  and  nonequilibrium  expansion 
process  are  evaluated  for  a  typical  nozzle  configuration  in  a  later 
subsection. 
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Nozzle  contours  for  inviscid  isentropic  flow  were  computed  using 
Rao's  thrust  optimization  technique  and  computing  machine  programs. 

The  details  of  the  procedures  are  given  in  Ref.  3.  These  nozzle  configu¬ 
rations  are  shown  in  Figures  5.  2  and  5.  3  for  the  two  nozzle  entrance 
conditions.  Superimposed  on  these  figures  are  lines  of  constant  nozzle 
impulse  coefficient  Cj^j .  The  representation  used  in  these  figures  allows 
us  to  choose  the  nozzle  configuration  that  yields  optimum  thrust  within 
a  prescribed  length. 

The  results  presented  in  Figures  5.  2  and  5.  3  pertain  to  only 
inviscid  flow.  The  influence  of  viscosity  is  assumed  to  affect  the  flow 
in  a  narrow  region  called  boundary  layer,  near  the  nozzle  wall.  Follow¬ 
ing  the  well  known  concepts  of  boundary  layer  theory  one  can  compute 
the  mass  flow  defect  and  the  momentum  defect  of  the  flow  confined  in 
the  boundary  layer  along  the  wall.  The  boundary  layer  is  assumed  com¬ 
pressible  and  turbulent.  The  boundary  layer  at  the  nozzle  entrance  section 
is  assumed  equal  in  thickness  to  the  boundary  layer  seven  feet  behind 
the  leading  edge  of  a  flat  plate.  Nozzle  wall  temperature  was  taken  into 
account  in  the  boundary  layer  equations  by  the  inclusion  of  a  parameter 


At  each  of  the  two  flight  conditions  ^  w  is  assigned  a  value  such  that 
the  nozzle  wall  temperature  would  be  2680°R.  The  growth  of  the  bound¬ 
ary  layer  along  the  nozzle  wall  is  computed  taking  into  account  the  local 
flow  conditions  obtained  from  the  inviscid  flow  computations.  The  details 

66 


CONFIDENTIAL 


FIGURE  5.  2 

NOZZLE  IMPULSE  COEFFICIENT  CHART  FOR  INVISCID  FLOW 
Mc  =  3.31,  Y  =  1.  105 
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FIGURE  5. 3 

NOZZLE  IMPULSE  COEFFICIENT  CHART  FOR  INVISCID  FLOW 
Mc  =  4.02,  Y  =  1.179 
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of  analysis  and  computational  procedures  for  the  boundary  layer  are 
given  in  Ref.  4.  Only  the  pertinent  results  are  briefly  discussed  here. 
From  the  boundary  layer  computations  the  displacement  thickness 
is  computed  along  each  nozzle  wall  contour  shown  in  Figures  5.  2  and 
5.  3  for  the  respective  flow  conditions  along  the  wall.  Following  the 
procedure  employed  in  wind  tunnel  nozzle  designs,  one  can  enlarge  the 
cross-sectional  areas  of  the  inviscid  nozzle  contour  by  the  amount  of 
displacement  thickness,  as  shown  in  Figure  5.4,  and  expect  to  have 
the  inviscid  flow  field  left  unaltered.  The  nozzle  contours  thus  corrected 
for  the  boundary  layer  displacement  thickness  are  shown  in  Figures 
5.  5  and  5.  6. 

From  the  inherent  assumption  of  constant  static  pressure  normal 
to  the  wall  within  the  boundary  layer,  the  wall  pressures  for  the  contours 
shown  in  Figures  5.  5  and  5.  6  can  be  obtained  from  wall  pressures  com¬ 
puted  for  inviscid  flow  in  nozzle  configurations  shown  in  Figures  5.  2 
and  5.  3.  Thus,  one  can  evaluate 


the  integral  of  pressure  forces  on  the  divergent  portion  of  the  nozzle  wall 
is  carried  out  from  up  to  the  exit  area  A^*,  which  now  includes  the 
boundary  layer  displacement  thickness  at  the  nozzle  exit.  The  frictional 
forces  on  the  nozzle  wall  contribute  to  a  decrement  in  thrust.  By  an 


*In  the  case  of  two-dimensional  nozzles,  cross-sectional  area  A  and 
cross-section  height  y  are  interchangeable  in  the  following  discussion. 
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FIGURE  5.4 

SKETCH  SHOWING  THE  ACCOMMODATION  FOR  BOUNDARY 
DISPLACEMENT  THICKNESS 
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appropriate  choice  of  friction  coefficient  applicable  to  the  exhaust 
nozzle  of  the  SCRJ  engine,  the  frictional  losses  can  be  evaluated.  A 
detailed  discussion  of  the  choice  of  the  friction  coefficient  is  given  in 
Ref.  3.  The  exit  impulse  of  the  nozzle  contours  shown  in  Figures  5.  5 
and  5.  6  can  now  be  computed  as 


Jfc  =  J5  + 


p 


Afc 

PdA  -  wall  friction  drag 


(5.8) 


The  nozzle  impulse  coefficient 

Cnj  =  ^  '  I 

can  represent  the  influence  of  nozzle  configuration  and  viscous  losses 
when  is  computed  according  to  the  above  definition.  The  values  of 
Cj^  thus  computed  for  the  respective  contours  at  the  two  flight  condi¬ 
tions  are  also  shown  in  Figures  5.  5  and  5.  6  by  drawing  constant 
lines. 


When  we  examine  nozzle  contours  of  various  lengths,  but  of  same 
exit  area  we  find  that  beyond  a  certain  length  the  decrement  in  nozzle 
impulse  coefficient  due  to  wall  friction  overcomes  the  gain  due  to  lower 
flow  divergence  losses.  Thus,  for  a  given  area  ratio  there  is  an  optimum 
value  of  attainable.  These  values  are  indicated  in  Figures  5.  5  and 
5.6  for  the  two  families  of  nozzle  contours  computed  for  the  two  flight 
conditions.  In  Figure  5.  7,  a  comparison  of  values  for  one- 
dimensional,  and  contoured  nozzles  is  shown.  The  influence  of  wall 
friction  losses  on  is  also  shown  in  the  figure. 
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5.  5  EFFECT  OF  NONUNIFORM  CONDITIONS  AT  NOZZLE  ENTRANCE 

The  impulse  of  the  exhaust  gases  at  the  nozzle  exit  would  depend 
upon  the  entrance  flow  conditions.  Holding  the  value  of  the  same,  the 
nonuniform  conditions  at  the  nozzle  entrance  can  be  examined  in  terms  of 
exit  impulse.  Thus,  the  effect  of  combustor  exit  conditions  are  accounted 
for  in  the  nozzle  performance  instead  of  the  combustor  performance.  Since 
the  nonuniformities  could  result  in  nonisentropic  and  nonisenthalpic  flows 
with  a  certain  amount  of  mixing,  any  rigorous  solution  of  the  problem  is 
not  feasible.  Also,  it  is  not  possible  to  consider  parametric  variation  of 
nozzle  configuration  with  parametric  variation  of  entrance  conditions. 
However,  examination  of  nozzle  flows,  with  a  few  suitably  prescribed 
entrance  conditions,  can  yield  an  indication  of  the  influence  of  nozzle  en¬ 
trance  nonuniformities. 

For  the  purpose  of  studies  reported  here  a  nozzle  configuration  of 
(L/Yg)  =  100  and  (Y^/Y,.)  =  25  is  chosen.  This  is  one  of  the  family  of  noz¬ 
zle  contours  designed  for  uniform  entrance  conditions  with  M,.  =  4.  02  and 
)(  =  1.  179  and  is  shown  in  Figure  5.8.  From  the  flow  analysis  by  the 
method  of  characteristics,  the  entire  flow  field  in  the  nozzle  is  known. 

In  order  to  apply  the  "stream-tube"  method  for  the  study  of  nonuniformity 
effects,  the  nozzle  entrance  section  is  divided  into  three  equal  stream 
tubes,  each  passing  equal  amounts  of  mass  flow.  Utilizing  the  flow  field 
computed  for  Y  =  1 .  179,  stream  lines  are  drawn,  as  shown  in  Figure 
5.8,  separating  the  nozzle  flow  into  three  stream  tubes,  referred  to  as 
a,  b,  and  c,  in  the  following  discussions.  The  exit  flow  of  the  nozzle  con¬ 
figuration  is  nonuniform.  For  example,  the  static  pressures  across 
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the  nozzle  exit  are  nondimensionalized  to  entrance  station  pressure  P,. 
and  exhibit  a  distribution  as  shown  in  Figure  5.  9.  The  exit  £low  of  the 
nozzle  is  not  axial  and  the  flow  directions  measured  from  the  axial  direc¬ 
tion  are  shown  in  Fig.  5.  10.  It  should  be  remembered  that  the  information 
presented  in  Figures  5.  8  to  5.  10  iB  obtained  from  field  computed  for 
*  1.  179  and  uniform  nozzle  entrance  Mach  number  M,.  =  4.  02. 

In  evaluating  nozzle  exit  flow  condition  for  nonuniform  entrance 
conditions  the  "stream -tube"  method  of  analysis  is  employed.  Two  typical 
cases  of  nonuniform  conditions  are  assumed  at  the  nozzle  entrance  and 
the  resulting  exit  flow  impulses  are  compared  with  the  case  of  uniform 
entrance  flow.  The  pressure  and  gas  velocity  at  the  nozzle  entrance 
for  all  three  cases  are  held  the  same  and  considered  uniform.  Fuel  mass 
fraction  is  varied  between  the  stream  tubes,  keeping  the  total  fuel  mass 
fraction  the  same.  The  various  flow  parameters  assumed  in  the  three 
stream  tubes  at  the  nozzle  entrance  are  described  in  Table  5.  1. 

In  evaluating  the  flow  conditions  at  the  nozzle  exit  the  following 
assumptions  are  made: 

1.  One -dimensional  flow  in  each  stream  tube  without  any  mixing  be¬ 
tween  adjacent  stream  tubes. 

9 

2.  The  expanding  gases  are  in  chemical  and  thermodynamic  equilibrium 
in  each  stream  tube. 

3.  The  cross-sectional  areas  of  the  stream  tubes  at  the  nozzle  exit 
projected  on  the  exit  plane  normal  to  the  axis,  as  shown  in  Figure 
5.  11,  sum  up  to  the  nozzle  exit  area. 

4.  The  variation  of  static  pressure  across  the  nozzle  exit  as  obtained 
from  stream  tube  analysis  is  similar  to  that  obtained  in  the  case  of 
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FIGURE  5.9 

STATIC  PRESSURE  VARIATION  AT  THE  NOZZLE  SHOWN 
IN  FIGURE  5.  8 
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FIGURE  5.10 

VELOCITY  DIRECTIONS  AT  THE  EXIT  OF  THE  NOZZLE 
SHOWN  IN  FIGURE  5.  8 


79 


APL  TDR  64-112 


CONFIDENTIAL 


TABLE  5.  1 


FLOW  CONDITIONS  ASSUMED  AT  THE  NOZZLE  ENTRANCE 


Stream- Tube 

Stream- Tube 

Stream- Tube 

a 

b 

c 

cross-section 

Hr-  —  0.  6667  to  1.  0 
*5 

X-  =  0.  3333  to.  6667 
Y5 

'«»—  =  0  to  0.  3333 
*5 

pressure 

0 .  444 /atm 

0.  444/  atm 

0.  444  atm 

0.  71 17  x  10°  cm/  sec, 

velocity 

0.  7117  x  10°  cm/sec 

0.  71 1 7  x  10°cm/  sec. 

Case  1 

temperature 

3148°K 

3148°K 

3148°K 

fuel  mas  8 
fraction 

0.  02833 

0.  02833 

0.  02833 

Case  2 

temperature 

3061°K 

3148°K 

3237°K 

fuel  mass 
fraction 

0. 042495 

0.  02833 

0.  014165 

Case  3 

temperature 

3237°K 

3148°K 

3061°K 

fuel  mass 
fraction 

0.014165 

0.  02833 

0. 042495 
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LATERAL  FORCE  |ROM  nth  STREAM  TUBE  ■ 

-  Pn  *n  ^1  «■  5^2.  )  SIN  9n  ♦  PhAn  TAN  0n  COS  0n 


FIGURE  5.  1 1 

SKETCH  OF  FLOW  CONDITIONS  IN  THE  N-TH  STREAM  TUBE  AT 
NOZZLE  EXIT 
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y  «  1.  179  and  shown  in  Figure  5.9,  i.e. ,  even  though  (P^/Pg) 
values  are  not  the  same  as  shown  in  Figure  5.9,  the  variation  of 
P^,  as  denoted  by  (P^/P^)  is  the  same.  The  subscript  w  indicates 
the  value  of  P^  at  the  nozzle  wall. 

5.  The  flow  directions  of  the  stream  tubes  at  the  nozzle  exit  are 
the  same  as  those  obtained  in  the  flow  field  computed  for 
■  1.  179  and  shown  in  Figure  5.  10. 

The  evaluation  of  flow  conditions  in  the  stream  tubes  at  the  nozzle 

exit  is  carried  out  by  an  iterative  process  because  of  the  stipulation  that 

the  total  cross-sectional  area  at  exit  is  equal  to  Y^.  After  determining 

the  pressures,  densities,  and  gas  velocities  in  the  respective  stream  tubes 

one  can  evaluate  the  exit  impulse  function  as 

L  ®  -t-  R  ,  {~(rnn\4n  +  Pfcn  Ac, n)  Cos 06 n + 

(5.9) 


Fin  A^Tart  06 -  wall  friction  drag^j 

The  subscript  n  refers  to  the  n1*1  stream  tube,  subscript  6  refers  to  the 
nozzle  exit  location,  and  9  denotes  the  direction  of  stream  tube  measured 
from  the  nozzle  axis,  ih  is  the  total  mass  flow  through  the  nozzle.  Ar  is 
the  cross-sectional  area  of  the  n*^1  stream  tube,  which  in  the  two-dimensional 
case  can  be  replaced  by  Y  .  The  condition  of  prescribed  value  of  exit 
area  is  satisfied  by  requiring 


aVc  CoS  6<»n 


-  A6 


(5.10) 


The  wall  friction  drag  in  the  above  Eq.  (5.9)  is  assumed  to  be  at  the  same 
value  as  computed  for  nozzle  flow  field  obtained  for  Y  -  1.  179. 
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In  a  similar  manner  one  can  determine  the  value  of  which  is 
simpler  because  the  flow  is  axial  and  the  nonuniformities  are  only  in  gas 
density.  The  nozzle  impulse  coefficient  C^j  computed  for  the  three  cases 
under  consideration  are  given  in  Table  5.  2.  It  should  be  noted  that  the 
Cj^  values  shown  parametrically  in  Figure  5.  6  are  the  results  of  compu¬ 
tations  for  constant  =  1.  179,  whereas  the  value  shown  for  Case  1  in  the 
Table  5.  2  is  for  equilibrium  flow.  For  the  sake  of  comparison  the  value 
of  obtained  for  this  nozzle  from  constant  X  computations  is  also 
shown  in  the  Table.  The  influence  of  nonuniformity  at  nozzle  entrance 
on  C-t.  is  indicated  in  Table  5.  2.  To  eliminate  the  direct  influence  of 
nozzle  area  expansion  ratio,  the  values  of  are  nondimensionalized  to 
j^.,  the  one -dimensional  value  for  the  area  ratio  25.  For  the  cases  1,2,. 
and  3  the  value  of  j^.  is  taken  from  one -dimensional  equilibrium  flow 
computation  corresponding  to  uniform  entrance  conditions  as  in  Case  1, 
and  the  ratios  obtained  are  shown  in  Table  5.  2.  In  comparing  the 

results  from  constant  ^  calculation  the  value  used  for  j^.  was  taken  from 
one -dimensional  flow  at  Y  =  1.  179.  The  corresponding  values  of 
^Nj^C*Nji^  are  also  shown  in  the  table. 

The  results  presented  in  Table  5.  2  indicate  that  the  loss  in  thrust 
due  to  the  nonuniform  entrance  conditions  of  the  Case  2  or  Case  3  is  the 
same.  When  we  consider  uniform  entrance  flow  with  equilibrium  expan¬ 
sion,  the  influence  of  the  particular  nozzle  contour  and  friction  drag  on 
the  wall  together  amounted  to  a  loss  of  0.47  percent  of  one -dimensional 
exit  impulse.  On  the  other  hand,  similar  evaluation  based  on  -  1.  179 
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TABLE  5.  2 

EFFECT  OF  NONUNIFORM  CONDITIONS  AT  ENTRANCE  OF 
NOZZLE  CONFIGURATION  OF  FIGURE  5.  8 


CNj 

CNj/CNji 

CEj 

Uniform  entrance 
conditions  -  Case  1 

0.05752 

0. 92017 

0. 99530 

Nonuniform  entrance 
conditions  -  Case  2 

0.05697 

0. 91137 

0. 99479 

Nonuniform  entrance 
conditions  -  Case  3 

0.  05697 

0. 91137 

0. 99479 

Uniform  entrance 
condition  --one- 
dimensional  nozzle 

0. 06251 

1.0 

1.0 

Uniform  entrance 
conditions --If*  1.  179 
(from  Figure  5.  6) 

0.  1033 

0.91424 

0. 99129 

Uniform  entrance 
conditions  -  -  one  - 
dimensional  nozzle, 
r=  1.  179 

(from  Figure  5.  7) 

0. 11299 

1.0 

1.0 
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computations  indicate  a  loss  of  0.  871  percent  of  one -dimensional  exit 
impulse  for  ^  =  1.  179.  The  influence  of  the  gas  properties  show  up 
more  pronounced  on  the  net  nozzle  thrust  as  indicated  by  values. 
For  constant  ^  =  1.  179  computations  a  value  of  0.  1033  is  obtained  of 
C^j.,  compared  with  only  0.  05752  when  we  consider  equilibrium  flow. 

5.  6  LATERAL  FORCES  DUE  TO  NONUNIFORM  ENTRANCE 
CONDITIONS 


Depending  upon  the  configuration  of  the  combustion  chamber  and 
fuel  injector  the  flow  conditions  at  nozzle  entrance  could  be  nonuniform. 
These  nonuniformities  at  nozzle  entrance  would  affect  the  downstream 
flow  field  in  the  nozzle,  giving  us  differences  in  stream  thrust  at  the 
nozzle  exit  and  consequently,  different  values  of  lateral  loads.  Variations 
in  this  lateral  load  will  reflect  as  changes  in  vehicle  total  lift.  For  the 
purpose  of  this  investigation  the  nozzle  configuration  and  entrance  condi¬ 
tions  used  in  Section  5.  5  are  considered.  Using  the  "stream  tube" 
method  of  analysis  and  the  assumptions  listed  in  Section  5.  5  the  flow 
conditions  at  the  nozzle  exit  section  are  determined.  Even  though  the  flow 
directions  of  the  individual  stream  tubes  at  the  nozzle  exit  are  considered 
unaffected  by  the  nozzle  entrance  nonuniformities,  the  changes  in  exit 
stream  thrust  of  individual  stream  tubes  yield  changes  in  the  lateral 
force  on  the  nozzle.  After  evaluating  the  appropriate  exit  flow  conditions 
the  lateral  force  can  be  computed  as 

-4-  ffl^nVfen  +  RnA4*i)Sw94„+  Pfen 

m  (V  '  /e  , 


(5.  11) 
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where  the  notation  used  is  the  same  as  that  in  Eq.  (5.  9).  The  lateral 
forces  computed  for  nozzle  entrance  conditions,  defined  as  Cases  1,  2, 
and  3  in  Table  5.  1  are  shown  in  Table  5.  3  after  normalizing  to  j^.. 
Although  the  lateral  force  is  only  a  small  fraction  of  exit  stream  thrust 
it  is  of  the  same  order  of  magnitude  as  the  nozzle  net  thrust,  i.e.  ,  by 
comparing  to  C^j  values.  Hence,  the  changes  in  the  lateral  loads  in 
Cases  2  and  3  from  that  in  Case  1  are  of  significant  value.  If  the  nozzle 
entrance  conditions  are  unsteady  in  nature,  one  can  expect  the  lateral 
force  on  the  nozzle  also  to  be  unsteady.  Even  though  periodic  variations 
in  the  net  thrust  of  an  SCRJ  engine  can  be  tolerated  to  some  extent, 
periodic  changes  in  the  lateral  force  can  be  dangerous,  affecting  the 
flight  dynamics  and  structural  integrity  of  the  engine.  Solution  of  un¬ 
steady  flow  field  in  a  two-dimensional  nozzle  is  not  feasible.  Examina¬ 
tion  of  results  presented  in  Table  5.  3  indicates  that  if  the  nozzle  entrance 
conditions  change  periodically  from  Case  2  to  that  of  Case  3,  there  could 
be  periodic  variations  in  lateral  force  with  an  amplitude  of  .  004  times 
j^..  Comparing  with  values  of  C^,  we  find  that  the  amplitude  of  lateral 
force  variation  can  be  5.  1  percent  of  net  nozzle  thrust,  whereas  the 
amplitude  of  variation  of  net  nozzle  thrust  is  only  0.  96  percent  (from 
Cases  2  and  3  in  Table  5.  3). 
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TABLE  5.  3 


LATERAL  FORCES  DUE  TO  NONUNIFORM  ENTRANCE  CONDITION 


Lateral  Force* 

j6i 


Uniform  entrance, 
equilibrium  flow-- 

Case  1  0.05752  -0.0673 

Nonuniform  entrance, 
equilibrium  flow- - 

Case  2  0.05697  -0.0654 

Nonuniform  entrance, 
equilibrium  flow-  - 

Case  3  0.05697  -0.0694 


*  The  computed  lateral  force  is  in  the  negative  Y  direction. 
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5.  7  EFFECT  OF  NONEQUILIBRIUM  FLOW 

As  the  exhaust  gases  in  an  SCRJ  engine  expand  through  the  exit 
nozzle,  a  certain  amount  of  chemical  recombination  occurs  among  the 
constituent  species  of  the  gases.  The  amount  of  recombination  that 
occurs  depends  upon  the  time -rate  of  the  expansion  process  and  the 
chemical  kinetic  rates  of  the  various  reactions  that  can  take  place  among 
the  various  species  in  the  exhaust  gases.  Because  of  the  energy  released 
during  the  chemical  recombination,  the  exit  momentum  of  the  exhaust 
gases  increases  as  a  function  of  the  amount  of  recombination  that  can 
occur.  On  the  other  hand,  for  a  given  nozzle  configuration  of  prescribed 
area  ratio  and  prescribed  inlet  flow  conditions,  the  amount  of  recombina¬ 
tion  increases  with  the  nozzle  size,  i.e.  ,  a  linear  dimension  of  the 
nozzle.  The  "sudden-freezing"  concept  as  proposed  by  Bray  in  Ref.  5 
and  extended  to  the  hydrogen  air  case  by  Kushida  (Ref.  6)  can  easily  be 
applied  to  evaluate  the  expansion  process  in  a  nozzle,  hence  the  exit 
impulse.  The  validity  of  the  "sudden -free zing"  point  concept  in  the 
exhaust  nozzle  performance  analyses  has  been  demonstrated  by 
Westenberg  (Ref.  7). 

With  complex  chemical  systems  in  exhaust  nozzle  flow  there 
are  several,  rather  than  one,  possible  freezing  points.  In  the  method 
of  Ref.  6,  the  freezing  of  the  concentrations  of  the  hydrogen  atom  H, 
the  oxygen  atom  O,  and  the  hydroxyl  radical  OH  is  examined.  The 
freezing  point  furthest  downstream  in  the  stream  tube  is  chosen  as  the 
freezing  point  of  reaction.  The  freezing  point  functions  are  expressed 
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in  the  form 
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where  k^,  k^>  k^>  and  k^  are  the  reaction  rate  coefficients  listed  in 
Table  5.  4  for  the  four  important  recombination  reactions  in  hydrogen- 
oxygen  combustion.  These  functions  are  evaluated  as  part  of  an  equi¬ 
librium  expansion  program  and  are  independent  of  the  nozzle  size.  At 
a  given  pressure  condition  the  freezing  point  function  with  the  largest 
absolute  value  is  chosen  as  the  critical  one.  In  the  present  problem, 
the  OH  radical  proved  to  be  the  most  important.  The  values  of 
(-dflnP/ds)OH  are  shown  on  a  log-log  plot  as  a  function  of  pressure  in 
Figure  5.  12  for  various  values  of  fuel  mass  fractions.  The  fuel  mass 
fractions  are  normalized  to  the  stoichiometric  fuel  mass  fraction  and 
are  shown  as  equivalence  ratios  in  Figure  5.  12. 

TABLE  5.  4 

REACTION  RATES  USED  IN  FREEZING  POINT  COMPUTATION 
Reaction 


3  2 

Rate  Coefficient  (cm  /mole)  /  sec 


H 

+ 

OH  -^►H20 

H 

+ 

H  -£4*»H2 

O 

+ 

(M 

O 

A 

o 

H 

+ 

o  -^-OH 

k2  =  5.4  x  1017/T 
18, 


k2  =  7.  5  x  10AO/T 
D1' 

15 


k 3  =  0.  67 x  1014(4500/T)1,  5 


k4=  1,0x10 
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PRESSURE  (atm) 


FIGURE  5;  12 

CRITERION  FOR  EQUILIBRIUM  IN  THE  NOZZLE  FLOW 
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At  any  particular  value  of  pressure  P,  if  the  expansion  along  a 

j  J?nP 

stream  line  in  a  nozzle  produces  a  higher  value  of-  M  y- ■  than  indicated 
in  Figure  5.  12,  the  chemical  reactions  are  considered  frozen  at  that 
location.  Subsequent  expansion  of  the  exhaust  gases  as  they  travel  along 
the  stream  line  can  be  computed  with  frozen  chemical  composition  ob¬ 
tained  at  the  "sudden-freezing"  point. 

For  the  investigations  reported  here,  the  nozzle  configuration 
shown  in  Figure  5.  8  is  considered.  From  the  known  flow  field  previously 
computed  for  *  4.  02  and  )(  *  1.  179,  two  stream  lines  are  shown  in 
Figure  5.  8.  Including  the  nozzle  wall  and  center  line  we  have  four 
stream  lines  along  which  the  pressure  variations  are  known  from 
Y  =  1.  179  computation  and  are  described  in  Figures  5.  13a  and  5.  13b. 

For  a  particular  size  of  the  nozzle  entrance,  say  =  1  ft.  ,  the  values 
of  -  can  be  computed  along  these  stream  lines.  These  values  for 

each  stream  line  are  shewn  in  Figure  5.  14,  against  pressure  on  a 
log -log  plot.  Comparing  the  conditions  attained  in  the  nozzle,  as  shown 
in  Figure  5.  14  with  the  "sudden -freezing"  criterion,  as  shown  in  Figure 
5.  12,  one  can  locate  along  each  stream  line  the  point  downstream  of 
which  the  expansion  process  can  be  assumed  to  occur  with  frozen  compo¬ 
sition,  as  shown  in  Figure  5.  15.  In  evaluating  the  impulse  at  the  nozzle 
exit  the  following  assumptions  are  made: 


1.  The  nozzle  flow  is  divided  into  three  stream  tubes  with  no 
mixing  occurring  across  the  stream -tube  walls. 

2.  One -dimensional  flow  is  assumed  in  each  stream  tube,  i.e.  , 
flow  conditions  are  uniform  across  each  stream  tube  even 
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P  (otm) 

FIGURE  5.  14 

GRADIENT  OF  PRESSURE  ALONG  EACH  STREAM  LINE  IN 
NOZZLE  SHOWN  IN  FIGURE  5.  8 
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though  there  could  be  variations  between  stream -tube  to 
stream -tube. 

3.  Equilibrium  expansion  is  assumed  in  each  stream-tube  up 
to  the  "sudden-freezing"  line  and  subsequent  expansion  is 
considered  to  be  with  frozen  composition.  This  location 
for  "sudden-freezing"  is  described  in  terms  of  an  average 
pressure  across  the  stream-tube  in  our  computations. 

4.  The  pressures  Pg  reached  by  the  nozzle  exit  in  each  indi¬ 
vidual  stream  tube  are  such  that  the  variation  of  the  pres¬ 
sure  Pg  across  the  exit  section  is  similar  to  that  obtained 
from  the  flow  field  computed  for  Y  *  1.  179  and  shown  in 
Figure  5.  9. 

5.  The  values  of  pressures  in  the  stream -tubes  are  adjusted 
such  that  the  sum  total  of  the  stream -tube  cross-sections 
projected  on  the  exit  plane,  as  shown  in  Figure.  5.  11  are 
equal  to  the  prescribed  exit  area  A^. 

6.  The  velocity  directions  of  the  flow  in  the  stream-tubes  are 
assumed  to  be  similar  to  the  flow  field  computed  for 

Y  =  1.  179  and  shown  in  Figure  5.  10. 

As  in  the  previous  Section  5.  5  the  evaluation  of  flow  condi¬ 
tions  in  the  stream-tubes  at  the  nozzle  exit  is  carried  out  by  an 
iterative  process  by  satisfying  Eq.  (5.  10)  so  that  prescribed  exit  area 
A^  is  obtained.  After  determining  the  pressures,  densities,  and  gas 
velocities  in  the  respective  stream -tube  the  exit  impulse  function  j^ 
is  computed  according  to  Eq.  (5.  9).  Because  the  conditions  at  the 
nozzle  entrance  are  assumed  uniform  the  value  of  jg  will  be  the  same 
as  that  computed  for  Case  1  in  Table  5.  2.  The  nozzle  impulse  coef¬ 
ficient  and  the  exit  impulse  coefficient  computed  for  this 
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nozzle  under  the  "sudden-freezing"  assumption  (referred  to  as  Case  4 
in  following  discussions)  are  shown  in  Table  5.  5  and  are  compared  with 
the  values  for  equilirbium  flow  conditions. 

TABLE  5.  5 

EFFECT  OF  INCOMPLETE  RECOMBINATION  IN 

NOZZLE  CONFIGURATION  OF  FIGURE  5.  8  -  Y  =  1  ft. 

5 


C.r 

Nj 

CNji 

CEj 

Equilibrium  flow 
(Case  1) 

0. 05752 

0.  92017 

0. 99530 

"Sudden-freezing" 
assumption  (Case  4) 

0.05648 

0.  90353 

0. 99432 

The  results  reported  here  indicate  that  the  losses  due  to  incomplete 
recombination  of  exhaust  gases  in  the  nozzle  configuration  and  size  con¬ 
sidered  can  be  estimated  to  be  0.  098  percent  of  j^.  by  utilizing  the  "sudden- 
freezing"  method.  Even  though  this  loss  is  small  in  terms  of  it  can 
be  a  significant  portion  of  net  nozzle  thrust,  as  can  be  observed  from 


5.  8  INTERACTION  OF  NOZZLE  LOSSES 

From  the  discussions  presented  in  the  previous  subsections  one 
can  enumerate  the  various  losses  that  can  occur  in  a  nozzle  of  prescribed 
area  ratio  and  length. 
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1.  Losses  due  to  nonuniform  exit  flow  (even  in  the  case  of  optimum 
thrust  contour  the  exit  flow  deviates  considerably  from  one- 
dimensional  value  for  the  area  ratio). 

2.  Losses  due  to  wall  friction  on  the  nozzle  wall  and«the  heat 
transfer  to  nozzle  wall. 

3.  Losses  due  to  nonuniform  conditions  at  the  entrance  to  the 
nozzle. 

4.  Losses  due  to  incomplete  recombination  of  the  exhaust  gases. 

The  first  item  in  the  above  is  independent  of  nozzle  size  and  to 
some  extent  depends  upon  the  properties  of  the  exhaust  gases,  such  as 
an  appropriate  value  of  )f  used  for  the  expansion  process.  The  latter 
three  items  can  be  evaluated  only  for  a  prescribed  nozzle  size  and  flow 
conditions  in  the  nozzle. 

The  above  itemized  losses  are  influenced  by  each  other  and  one 
can  carry  out  the  computations  by  suitable  assumptions  only  for  a  given 
nozzle  flow  field,  as  described  in  detail  in  Sections  5.  5  and  5.  7.  . 
However,  one  can  examine  whether  these  losses  can  be  computed  indi¬ 
vidually  so  that  tedious  computations  can  be  avoided.  The  following 
discussion  pertains  to  checking  the  validity  of  adding  up  individual  losses 
computed  from  quasi  one -dimensional  flows  to  obtain  the  combined  loss 
in  a  given  nozzle  configuration. 

The  first  of  the  various  losses  itemized  above  can  be  computed 
for  a  given  nozzle  configuration  (area  ratio  and  length)  based  on  optimum 
thrust  contour  and  an  appropriate  value  of  constant  V.  For  the  nozzle 
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configuration  described  in  Figure  5.8  this  flow  divergence  loss  is  8.  58 
percent  of  CNj.,  with  CN„  computed  for  ^=1.  179  and  (A^/Ag)  =  25. 

The  second  item,  i. e.  ,  loss  due  to  frictional  drag  on  the  nozzle 
wall  was  computed  as  1.  51  percent  of  C^,  based  on  wall  conditions 
obtained  from  constant  ^  =  1.  179  computations.  One  may  say  that 
minor  variations  of  nozzle  contour  and  equilibration  of  exhaust 
products  would  not  alter  the  frictional  drag  loss  appreciably,  provided 
the  linear  dimensions  are  kept  the  same.  The  influence  of  heat 
transfer  to  nozzle  wall  is  not  considered  in  the  present  evaluation. 

The  last  item,  i.  e.  ,  loss  due  to  incomplete  recombination  can 
be  easily  evaluated  for  quasi  one -dimensional  flow  at  the  nozzle  exit. 

The  exit  flow  of  the  nozzle  can  be  assumed  at  uniform  pressure  and 
axial  direction  in  place  of  preassigned  distributions  as  in  Figures  5.9 
and  5.  10.  The  flow  through  the  nozzle  is  divided  into  three  stream 
tubes  as  in  Section  5.  7  to  obtain  "sudden-freezing"  locations  in  the 
individual  stream  tubes.  However,  uniform  pressure  is  assumed  at 
the  nozzle  exit  and  the  condition  of  assigned  value  of  A^  is  imposed. 
Because  the  amount  of  recombination  occurring  in  each  stream  tube  is 
different,  the  exit  specific  impulse  is  not  uniform  across  the  nozzle  exit. 
The  result  of  such  computations  shows  that  for  the  assumption  of 
uniform  pressure  at  the  exit  plane,  the  loss  indicated  for  incomplete 
recombination  is  1.  55  percent  of  Cj^. 
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TABLE  5.6 

COMPARISON  OF  LOSSES  ESTIMATED  BY 
QUASI  ONE -DIMENSIONAL  ANALYSIS  AND 
MORE  EXACT  ANALYSIS 

(The  losses  denoted  as  percent  of  Cj^) 


Case  2  Case  4 

Contour  loss  8.  58  8. 58 

Friction  loss  1.51  1.51 

Nonuniform  entrance  loss  0. 66  0 

Incomplete  recombination 

loss  0  1.55 

Estimated  loss  10.75  11.64 

(Quasi  one -dimensional) 

Exact  analysis  8. 86  9. 65 

(from  Tables  5.  2  &  5.  5) 
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The  third  item,  i.  e.  ,  losses  due  to  nonuniformities  at  nozzle 
entrance  can  also  be  evaluated  in  terms  of  one -dimensional  exit  flow. 
The  approach  taken  is  similar  to  the  one  described  above  except  that 
equilibrium  flow  is  assumed  in  each  stream  tube.  The  results  of  such 
computations,  for  the  two  types  of  nonuniform  entrance  conditions 
described  as  Cases  2  and  3  of  Table  5.  1  show  that  for  the  assumption 
of  uniform  exit  pressure  the  loss  due  to  nonuniform  entrance  conditions 
is  0.  66  percent  of  Cj^. 

The  various  losses  computed  from  quasi  one -dimensional 
analysis  are  shown  in  Table  5.  6.  Let  us  consider  the  summation  of 
the  losses.  For  the  nozzle  flow  conditions  described  as  Case  2,  the 
quasi  one -dimensional  estimates  indicate  a  loss  of  10.  75  percent  of 
Cj^  compared  with  a  slightly  lower  value  of  8.  86  percent  obtained 
from  the  more  exact  analysis.  For  the  nozzle  flow  conditions  described 
as  Case  4  the  quasi  one -dimensional  estimates  indicate  a  loss  of 
1 1.  64  percent  of  Cj^  compared  with  9.  65  percent  obtained  from  the 
more  exact  analysis. 

From  the  numerical  computational  results  and  discussion  pre¬ 
sented  in  this  subsection,  one  can  conclude  that  even  though  the  losses 
in  a  nozzle  depend  upon  the  particular  flow  field,  one  can  make  con¬ 
servative  estimates  from  quasi  one -dimensional  computations. 

5.  9  LIFT  FORCE  DUE  TO  ASYMMETRIC  CONFIGURATION 

Asymmetric  nozzle  configurations  have  appeared  with  consid¬ 
erable  frequency  in  supersonic  combustion  ramjet  conceptual  studies. 
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This  is  due  partly  to  the  problems  involved  with  integrating  the  engine 
with  an  airframe  and  partly  to  the  asymmetry  of  low  loss  inlets.  A 
lateral  force  is  generated  by  an  asymmetric  nozzle  which,  if  the 
nozzle  is  properly  oriented,  will  contribute  to  the  vehicle  lift. 

One  particular  type  of  asymmetric  nozzle  was  studied  in  Ref.  8. 
This  was  the  one  wherein  the  lower  wall  of  a  symmetric  nozzle  is  cut 
back  shorter  than  the  upper.  For  this  type  asymmetry,  the  thrust 
and  lift  forces  may  be  treated  independently  as  long  as  the  last  left¬ 
running  characteristic  from  the  lower  wall  does  not  intersect  the 
upper  wall,  as  shown  in  Figure  5.  16. 

From  the  flow  fields  computed  for  constant  If  in  several 
uniform  exit  flow  nozzles,  the  ratios  of  lift  force  to  axial  component 
of  exit  stream  thrust,  obtained  for  various  amounts  of  lower  wall 
truncation,  are  shown  in  Figure  5.  17.  The  ambient  pressure  is 
assumed  zero,  giving  us  the  lift  force  directly  from  the  differential 
of  pressure  forces  acting  on  the  upper  and  lower  walls. 

To  indicate  the  loss  in  thrust  due  to  lower  wall  truncation, 
the  changes  in  the  exit  stream  thrust  are  shown  in  Figure  5.  18.  The 
specific  curves  (previously  described  in  Ref.  8)  obtained  for  each 
area  ratio  nozzle  crossed  one  another,  hence  have  been  shown  as  a 
band  for  simplicity  in  both  Figures  5.  17  and  5.  18.  It  should  be 
remembered  that  the  normalizing  factor  used  in  these  figures  is  the 
exit  stream  thrust,  whereas  the  net  engine  thrust  will  be  only  a 
fraction  of  this  due  to  the  inlet  momentum  drag.  An  examination  of 
Figures  5.  17  and  5.  18  indicates  that  considerable  amounts  of  lift 
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FIGURE  5.  16 

TRUNCATION  OF  LOWER  WALL  IN  A  TWO-DIMENSIONAL 
SYMMETRIC  EXIT  NOZZLE 
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FIGURE  5.17 

LATERAL  FORCE  DEVELOPED  BY  LOWER 
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LOSS  IN  EXIT  STREAM  THRUST  DUE  TO  LOWER 
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force  can  be  obtained  with  negligible  loss  in  thrust  performance. 

Consideration  of  nozzle  weight  and  envelope  restrictions  also  would 

favor  the  use  of  lower  wall  truncation. 

5.10  CONCLUSIONS 

Based  on  the  nozzle  flow  studies  and  the  results  presented 

here  the  following  conclusions  can  be  made. 

1.  Even  minor  gains  in  exit  stream  thrust  are  of  significant 
importance  in  the  case  of  an  SCRJ  engine  exit  nozzle,  hence 
it  is  profitable  to  use  nozzle  configurations  optimized  for 
thrust. 

2.  The  properties  of  exhaust  gases  employed  in  the  analysis  of 
nozzle  flow  (i.  e.  ,  constant  if ,  or  equilibrium  flow  assump¬ 
tions)  affects  the  nozzle  impulse  coefficient  considerably. 

3.  Frictional  drag  on  the  nozzle  wall  is  of  significant  value 
such  that  one  finds  an  optimum  length  of  a  given  nozzle 
area  expansion  ratio. 

4.  In  evaluating  the  losses  due  to  nonuniform  entrance  conditions 
or  incomplete  recombination  within  the  nozzle,  the  assumptions 
made  regarding  the  exit  flow  conditions  play  an  important 
role. 

5.  The  various  losses  in  a  nozzle  attributable  to  different  sources 
are  interdependent  and  cannot  be  summed  up  from  separate 
estimates.  The  sum  total  of  the  losses  has  to  be  evaluated 
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for  the  flow  field  occurring  in  the  nozzle  being  examined. 

6.  Asymmetric  nozzle  configurations  can  profitably  be  employed 
to  yield  large  amounts  of  lift  force  with  negligible  loss  in  nozzle 
thrust. 

7.  Nonuniform  conditions  at  nozzle  entrance  yield  lateral  forces 
on  the  nozzle,  which  vary  with  the  nature  of  the  nonuniformity. 
For  the  particular  types  of  nonuniformities  studied,  changes  in 
the  lateral  force  can  be  five  times  the  changes  in  the  net  nozzle 
thrust. 

5.  1 1  RECOMMENDATIONS 

The  exit  nozzle  studies  presented  in  this  section  and  the  limited 
number  of  computations  carried  out  showed  the  need  for  further  investi¬ 
gations  in  certain  areas.  Based  on  our  studies  the  following  items  can 
be  recommended  for  further  examination  to  provide  a  sound  basis  for 
nozzle  performance  estimates. 

1.  When  we  consider  nozzles  of  limited  length,  the  exit  flow  is 

necessarily  nonuniform  and  the  consequent  loss  in  thrust  is  a 
strong  function  of  the  exhaust  gas  properties  used  in  the  flow 
field  analysis.  One  is  tempted  to  use  an  appropriate  constant 
value  of  Y .  However,  the  particular  choice  of  Y  needs 
caution.  For  example,  the  flow  divergence  loss  obtained  for 
equilibrium  flow  does  not  apply  for  the  case  of  flow  with  in¬ 
complete  recombination.  Further  studies  on  the  appropriate 
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choice  of  %  for  evaluation  of  flow  field  is  recommended. 

2.  In  utilizing  "stream-tube"  method  for  evaluating  the  exit  stream 
thrust  of  a  nozzle  of  limited  length,  an  assumption  regarding  the 
nature  of  the  exit  flow  conditions  is  necessary.  However,  the 
results  appear  to  depend  upon  the  assumptions  made.  Further 
analytical  and  experimental  effort  is  recommended  to  guide  the 
investigator  in  the  choice  of  the  assumptions. 

3.  Periodic  nature  of  lateral  forces  caused  by  periodic  non¬ 
uniformities  at  nozzle  entrance,  require  further  study. 

4.  The  use  of  "sudden-freezing"  concept  coupled  with  "stream- 
tube"  method  appears  attractive  for  nozzle  flow  analysis  in  the 
case  of  incomplete  recombination  of  exhaust  gases.  However, 
the  validity  of  such  an  approach  remains  to  be  checked  against 
detailed  flow  field  computations  utilizing  the  method  of  charac¬ 
teristics,  incorporating  finite  rates  of  reactions  among  the 
species  of  the  exhaust  gases. 

5.  In  the  parametric  study  of  optimum  thrust  nozzle  configurations, 
it  is  not  necessary  to  confine  oneself  to  the  use  of  constant  Y  . 
Provided  the  "sudden-freezing"  concept  is  valid,  the  last 
characteristic  line  (in  most  cases)  would  lie  in  the  region  of 
flow  with  frozen  gas  composition.  Rao's  method  can  then  be 
applied  and  the  contour  optimization  can  be  extended  to  cover 
nozzles  where  recombination  is  incomplete.  Further  investi¬ 
gation  of  this  approach  is  recommended. 
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